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1 Introduction 

Homomorphisms between relational structures are not only fundamental mathematical objects, but are 
also of great importance in an applied computational context. Indeed, constraint satisfaction problems, 
a wide class of algorithmic problems that occur in many different areas of computer science such as 
artificial intelligence or database theory, may be viewed as asking for homomorphisms between two re- 
lational structures BFV98II . In a logical setting, homomorphisms may be viewed as witnesses for positive 
primitive formulas in a relational language. As we shall see below, homomorphisms, or more precisely 
the numbers of homomorphisms between two structures, are also related to a fundamental computational 
problem of statistical physics. Homomorphisms of graphs are generalizations of colorings, and for that 
reason a homomorphism from a graph G to a graph H is also called an H-coloring of G. Note that 
if H is the complete graph on /c-vertices, then an //-coloring of a graph G may be viewed as a proper 
/c-coloring of G in the usual graph theoretic sense that adjacent vertices are not allowed to get the same 
color. 

It is thus no surprise that the computational complexity of various algorithmic problems related to 
homomorphisms, in particular the decision problem of whether a homomorphism between two given 
structures exists and the counting problem of determining the number of such homomorphisms, have 
been intensely studied. (For the decision problem, see, for example, [B KN09I IBul06l IBKJ051 IGro07[ 
HN90]. References for the counting problem will be given in Section[3] Other related problems, such as 
optimization or enumeration problems, have been studied, for example, in BAus07[ IBDGM091 IDJKK081 
Rag08[[SS07).) 



In this article, we are concerned with the complexity of counting homomorphisms from a given 
structure A to a fixed structure B. Actually, we are mainly interested in a generalization of this problem 
to be introduced in the next section. We almost exclusively focus on graphs. The first part of the article, 
consisting of the following two sections, is a short survey of what is known about the problem. In the 
second part, consisting of the remaining Sections [4]|9j we give a proof of a theorem due to Bulatov and 
the first author of this paper [BG05 ], which classifies the complexity of partition functions described by 
matrices with non-negative entries. The proof we give here is essentially the same as the original one, 
with a few shortcuts due to IIThu09ll . but it is phrased in a different, more graph theoretical language that 
may make it more accessible to most readers. 



'supported in part by Marie Curie Intra-European Fellowship 271959 at the Centre de Recerca Matematica, Bellaterra, 
Spain 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



Qd — (D O^C® ® — © 

Figure 1: The graphs I, J, K% 

Ml) = (j J) A(J) = (J *) = (J J) A(A- S ) = (J J lj 

Figure 2: The adjacency matrices of the graphs I, K2, ks 

2 From Homomorphisms to Partition Functions 

For a fixed graph H we let Zh be the "homomorphism-counting function" that maps each graph G to 
the number of homomorphisms from G to H. Several well-known combinatorial graph invariants can 
be expressed as homomorphism counting functions, as the following examples illustrate: 

Example 2.1. Let I be the first graph displayed in Figure^ and let G be an arbitrary graph. Remember 
that an independent set (or stable set) of G is set of pairwise nonadjacent vertices of G. For every set 
S C V(G), we define a mapping hs : V{G) — > V(I) by /15(f) = 2 if v G 5 and hg{v) = 1 otherwise. 
Then hs is a homomorphism from G to I if and only if S is an independent set. Thus the number Zt(G) 
of homomorphisms from G to I is precisely the number of independent sets ofG. 

Example 2.2. For every positive integer k, let be the complete graph with vertex set [k] := {1, . . . , k} 
(see Figure^. Let G be a graph. Recall that a (proper) ^-coloring of G is a mapping h : V(G) — > [k] 
such that for all vw G E(G) it holds that h(v) 7^ h(w). Observe that a mapping h : V(G) —> [k] is 
a k-coloring of G if and only if it is a homomorphism from G to K\.. Hence Zx k (G) is the number of 
k-colorings of G. 

Unless mentioned otherwise, graphs in this article are undirected, and they may have loops and parallel 
edges. Graphs without loops and parallel edges are called simple. We always assume the edge set and 
the vertex set of a graph to be disjoint. The class of all graphs is denoted by Q. A graph invariant is 
a function defined on Q that is invariant under isomorphism. The adjacency matrix of a graph H is the 
square matrix A := A(H) with rows and columns indexed by vertices of H, where the entry A V)W at 
row v and column w is the number of edges from v to w. Figure [2] shows the adjacency matrices of the 
graphs in Figure [T] For all graphs G, H, we define a homomorphism from G — > H to be a mapping 
h : V(G) U E(G) -> V(H) U E(H) such that for all v G V(G) it holds that h(v) G V(H) and 
for all edges e G E(G) with endvertices v, w it holds that h(e) G E(H) is an edge with endvertices 
h(v), /i(u>)[j] The following observation expresses a homomorphism counting function Zh in terms of 
the adjacency matrix of H: 

Usually, homomorphisms from G to H are defined to be mappings g : V(G) — > V(H) that preserve adjacency. A 
mapping g : V(G) — > V(H) is a homomorphism in this sense if and only if it has an extension h : V(G) U E(G) — > 
V(H) U E(H) that is a homomorphism as defined above. Thus the two notions are closely related. However, if H has parallel 
edges, then there a different numbers of homomorphisms for the two notions. 
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Observation 2.3. Let H be a graph and A := A(H). Then for every graph G, 

z H {g)= yi n 4tw.«x(«)- w 

cr:V(G)^V(H) ee_B(G) with 
endvertices v,w 

To simplify the notation, we write j~J instead of in similar expressions in the following. 

vw£E(G) eSE(G) with 

endvertices v,w 

By convention, the empty sum evaluates to and the empty product evaluates to 1. Thus for the empty 
graph we have Z H {%) = 1 for all H and Z $ (G) = for all G / 0. 

Example 2.4. Consider the second graph J displayed in Figure [7] and let G be an arbitrary graph. 
Zj(G) is a weighted sum over all independent sets of G: For all sets S,T C V(G) we let e(S, T) be 
the number of edges between S and T. Then 

Zj(G)= 2 e{S ' V{G)XS) - 

SCV(G) 
independent set 

Equation ( |2.1| ) immediately suggests the following generalization of the homomorphism counting 
functions Zh'- For every symmetric n x n matrix A with entries from some ring § we let Za ■ Q — > § 
be the function that associates the following element of S with each graph G = (V, E): 

z A (G) = U Mvu^y w 

a:V— >[n] vw&E 

We call functions Za, where A is a symmetric matrix over S, partition functions over S. (All rings in 
this paper are commutative with a unit. S always denotes a ring.) 

With each n x n matrix A we associate a simple graph H(A) with vertex set [n] and edge set 
{ij | Aij 7^ 0}. We may view A as assigning nonzero weights to the edges of G(A), and we may 
view Za(G) as a weighted sum of homomorphisms from G to H(A), where the weight of a mapping 
a : V — > [n] is 

ua(G, ct) := Y\ At(i>),<7(u>)- 

vw£E 

Homomorphisms are precisely the mappings with nonzero weight. Inspired by applications in statistical 



physics (see Section 2.1 ), we often call the elements of the index set of a matrix, usually [n], spins, and 



we call mappings a : V — )■ [n] assigning a spin to each vertex of a graph configurations. 

Example 2.5. Recall that a graph G is Eulerian if there is a closed walk in G that traverses each edge 
exactly once. It is a well-known theorem, which goes back to Euler, that a graph is Eulerian if and only 
if it is connected and every vertex has even degree. Consider the matrix 



U 



1 -1 
-1 1 



It is not hard to show that for every N-vertex graph G we have Zjj(G) = 2 N if every vertex of G has 
even degree and Zjj(G) = otherwise. Hence on connected graphs, [1/2 N ) ■ Zjj is the characteristic 
function of Euler ianicity. 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



Example 2.6. Consider the matrix 

B :i 1 



1 -1, 

Let G be a graph. Then for every a : V(G) — > [2] it holds that 

1 if the induced subgraph G [a -1 (2)] has an even number of edges, 
— 1 otherwise. 



u B (G,a) -- 
It follows that for every N -vertex graph G, 



1 -Z B (G) + 2 N -i 

is the number of induced subgraphs of G with an even number of edges. 

Example 2.7. Recall that a cut of a graph is a partition of its vertex set into two parts, and the weight 
of a cut is the number of edges from one part to the other. A maximum cut is a cut of maximum weight. 
Consider the matrix 

1 X N 
.X 1 



C :-- 



over the polynomial ring 7h\X\. It is not hard to see that for every graph G, the degree of the polynomial 
Zc(G) is the weight of a maximum cut of G and the leading coefficient the number of maximum cuts. 

Graph polynomials present another important way to uniformly describe families of graph invariants. 
Examples of graph polynomials are the chromatic polynomial and the. flow polynomial. Both of these 
are subsumed by the bivariate Tutte polynomial, arguably the most important graph polynomial. The 
following example exhibits a relation between the Tutte polynomial and partition functions. 

Example 2.8. Let G = (V, E) be a graph with N vertices, M edges, and Q connected components. 
For a subset F C E, by q(F) we denote the number of connected components of the graph (V, F). The 
Tutte polynomial of G is the bivariate polynomial T(G; X, Y) defined by 

T(G;X,Y) = £ (X - l) q(F) - Q ■ (Y - 1)\ F \~ N +^ F ). 

FCE 

It is characterized by the following contraction-deletion equalities. For an edge e £ E, we let G \ e be 
the graph obtained from G by deleting e, and we let G/e be the graph obtained from G by contracting 
e. A bridge of G is an edge e £ E such that G\e has more connected components than G. A loop is an 
edge that is only incident to one vertex. 



'l ifE(G) 



T(G;X,Y) = { 



X ■ T(G \ e; X, Y) if e G E(G) is a bridge, 

Y ■ T(G/e; X, Y) if e G E(G) is a loop, 

T(G \ e; X, Y) + T(G/e; X, Y) if e G E(G)is neither a loop nor a bridge. 



Let r, s G C. It can be shown that the partition function of the n x n matrix A{n, r, s) with diagonal 
entries r and off-diagonal entries s satisfies similar contraction-deletion equalities, and this implies that 
it can be expressed in terms of the Tutte polynomial as follows: 

Z A{n ,r, s) {G) = s M - N +Q ■ (r - -rfl-T (g; * + ^ >~ ^ ^ . (2.3) 



2.1 Partition functions in statistical physics 
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This implies that for all x, y G C such that n := (x — 1) • (y — 1) is a positive integer, it holds that 

T (G; x, !/) = (!/- 1)°"^ • • Z A{n>yA) (G). 

Example 2.9. For simplicity, in this example we assume that G = (V, E) is a simple graph, that is, a 
graph without loops and without multiple edges. For every positive integer k, a fc-flow in G is a mapping 
f : V — >■ Zfe (the group of integers modulo k) such that the following three conditions are satisfied: 

(i) f(v, w) = Ofor all v,w G V with vw E; 

(ii) f(v,w) = —f(w,v)forallv,w G V; 
(Hi) ^2 w ) = °/° r °U v <EV. 

wdV with 

The k-flow f is nowhere zero if f(v, w) ^ Ofor all vw G E. Let F(G, k) be the number of nowhere 
zero k-flows of G It can be shown that 

F(G, k) = (-i)M-N+Q . t(G; 0, 1 - k) = k~ N ■ Z A{ktk _ lt _ x) {G) , 

where T denotes the Tutte polynomial and A(k, k — 1, 1) the k x k matrix with diagonal entries (k — 1) 
and off-diagonal entries — 1. 

2.1 Partition functions in statistical physics 

The term "partition function" indicates the fact that the functions we consider here do also have an 
origin in statistical physics. A major aim of this branch of physics is the prediction of phase transitions 
in dynamical systems from knowing only the interactions of their microscopic components. In this 
context, partition functions are the central quantities allowing for such a prediction. As a matter of fact 
many of these partition functions can be described in the framework we defined above. 

Let us see an example for this connection — the partition function of the Ising model. Originally 
introduced by Ising in 1925 [Isi25] this model was developed to describe the phase transitions in fer- 
romagnets. For some given graph G, the model associates with each vertex v a spin a v which may be 
either +1 or —1. Then the energy of a state a is given by the Hamiltonian defined by 

H(a) =~JJ2 u ^ (2 " 4) 

uv£E 

where —Ja u a v is the contribution of the energy of each pair of nearest neighbor particles. Let T denote 
the temperature of the system and k be Boltzmann's constant, define j3 = (kT)^ 1 . Then, for a graph 
G = (V, E) with N vertices, M edges, and Q connected components, we have 

Z(G,T) = ]T e-^W 

cr:V-*{+l,-l} 

We straightforwardly get Z(G, T) = e l3JM Z A (G) for the matrix 



A = A(T) = 




(2.5) 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



An extension of this model to systems with more than two spins is the n-state Potts model, whose 
partition function satisfies 



J Potts I 



In fact, for n = 2 and v = e 2/3J - 1 we see that Z(G, T) = e@ JM Z PottSi (G; n, v). Moreover, this model 
can be seen as a specialization of the Tutte Polynomial. Expanding the above sum over connected 
components of G, we obtain 



Z Potts (G;n,v) = y~] 



n l(A) v \A\ 



ACE 

whence it is not difficult to see that, if (X - 1)(Y - 1) G N, 

T(G; X, Y) = (X- 1)-^ E \Y - 1)~ N Z Potts (G; (X - 1)(Y — 1),Y — 1). 



Since Z Potts (G; n, v) = ZM n ,v+i,\){G) trns relation is actually a special case of equation (23). 



2.2 Which functions are partition functions? 

In this section, we shall state (and partially prove) precise algebraic characterizations of partition func- 
tions over the real and complex numbers. 

A graph invariant / : Q ->■ S is multiplicative if /(0) = 1 and f(G ■ H) = f(G) ■ f(H). Here 
denotes the empty graph and G-H denotes the disjoint union of the graphs G and H. An easy calculation 
shows: 



Observation 2.10. All partition functions are multiplicative. 



To characterize the class of partition functions over the reals, let / : Q — > R be a graph invariant. 
Consider the(infinite) real matrix M = (Mg,h)g,h&g with entries Mq^h '■= f(G • H). It follows from 
the multiplicativity that if / is a partition function then M has rank 1 and is positive semidennite. (Here 
an infinite matrix is positive semidefinite if each finite principal submatrix is positive semidennite.) The 
criterion for a graph invariant being a partition function is a generalization of this simple criterion. Let 
k > 0. A k-labeled graph is a graph with k distinguished vertices. Formally, a fc-labeled graph is a pair 
(G,4>), where G G Q and <fi : [k] — > V(G) is injective. The class of all /c-labeled graphs is denoted 
by Ok- F° r two fc-labeled graphs (G, 4>), (H, ijj) G Q^, we let (G, 4>) ■ (H, ip) be the /c-labeled graph 
obtained from the disjoint union of G and H by identifying the labeled vertices <f>(i) and ip(i) for all 
i G [k] and keeping the labels where they are. We extend / to Qk by letting f(G, 4>) := f{G) and define 
a matrix 

by letting M(f, k)( G ^ ih,^) = f{{G, <j>) • (H, -0)). Note that if we identify 0-labeled graphs with plain 
graphs, then M(f, 0) is just the matrix M defined above. The matrices M(f, k) for k > are called 
the connection matrices of /. Connection matrices were first introduced by Freedman, Lovasz, and 
Schrijver [FLS07] to prove a theorem similar to the following one (Theorem [2T4] below). 



Theorem 2.11 (Schrijver |Sch09|). Let f : Q -^-Mbe a graph invariant. Then f is a partition function 
if and only if it is multiplicative and all its connection matrices are positive semidefinite. 



2.2 Which functions are partition functions ? 
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We first give the simple proof of the forward direction. A proof sketch for the backward direction will 
be given at the end of this subsection. Suppose that / = Za for a symmetric matrix A G M nxn . For all 
mappings \ '■= [k] — > [n] and A;-labeled graphs (G, 4>) G Gk we let 

Za, x (G, (j)) := ^2 ] J At(v),<t(to)- 

<x:V(G)->[n] vweE{G) 
a(4>{i))=x{i) foraUie[fc] 

Note that for (H, ip) G <5 fc we have 

^, X ((G,0) • = ^,x(G,0) • Z A)X {H^). 

Hence the matrix M(Z A ,xi k) with entries M(Z A ,xi k)(G,c/>),(H,i>) = Za,x{{G, 4>) ■ {H, ifi)) is positive 
semidefinite. Furthermore, 

f{G,<j>) = Z A (G) = Z Ax(G,<l>), 

and thus Mf^ is the sum of n k positive semidefinite matrices, which implies that it is positive semidef- 
inite. Incidentally, the same argument shows that the row rank of the kth connection matrix M(Za, k) 
of a partition function of a matrix A G M. nxn is at most n k . 

Schrijver also obtained a characterization of the class of partition functions over the complex num- 
bers, which looks surprisingly different from the one for the real numbers. In the following, let S = C or 
C[X] for some tuple X of variables. For a symmetric n x n matrix A G § nxn we define the "injective" 
partition function Y A : Q — > E> by 

Y A (G) := ^ Y[ A t{v),t(w)- 

r:V(G)->[n] vweE(G) 
injective 

Note that Y A (G) = for all G with \V{G)\ > n. For a graph G and a partition P of V(G), we let G/P 
be the graph whose vertex set consists of the classes of P and whose edge set contains an edge between 
the class of v and the class of w for every edge vw G E(G). Note that in general G/P will have many 
loops and multiple edges. For example, if P has just one class, then G/P will consist of a single vertex 
with |-E(G)| many loops. We denote the set of all partitions of a set V by n(V), and for a graph G we 
let 11(G) := n(V(G)). For P, Q G U(V), we write P < Q if P refines Q. Then we have 

Z A (G)= J2 y a(G/P). (2.6) 
Pen(G) 

We can also express Y A in terms of Za- For every finite set V there is a unique function \i : IT (V) — > Z 
satisfying the following equation for all P G n(V): 

E m(Q) = {; if th P = rV ' (2 " 7) 
q^) 1° otherwise - 

Q<P 

Here Ty denotes the trivial partition {{v} \ v G V}. The function /i is a restricted version of the Mobius 
function of the partially ordered set n(V) (for background, see for example [Aig07]). Now it is easy to 
see that 

Y A (G)= E KP)-Za(G/P). (2.8) 
Pen(G) 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



We close our short digression on Mobius inversion by noting that for all sets V with |V| =: k we have 
the following polynomial identity: 

Y fi(P)-X^ =X(X -1)---(X -k + l). (2.9) 

Pen(v) 

To see this, note that it suffices to prove it for all X G N. So let XeN, and let A be the (X x X)- 
identity matrix. Furthermore, let G be the graph with vertex set V an no edges. Then Ya = X(X — 
1) • • • (X - k + 1) and Z A (G/P) = X\ p \ for every partition P G U(V), and (|Z9]) follows from Q. 



Now we are ready to state Schrijver's characterization of the partition functions over the complex 
numbers. 

Theorem 2.12 (Schrijver |Sch09|). Let f : Q — > C be a graph invariant. Then f is a partition function 
if and only if it is multiplicative and 

Y KP)-f(G/P) = (2.10) 
Pen(G) 

for all G G Q with \V(G)\ >\f(Ki)\. 

To understand the condition |V(G)| > \f(Ki)\, remember that K\ denotes the graph with one vertex 
and no edges and note that for every n x n matrix A it holds that Za(K{) = n. 



Proof of Theorem 2.12 (sketch). The forward direction is almost trivial: Suppose that / = Za for a 
symmetric matrix A G C nxn . Then / is multiplicative by Observation 2.10 and we have 



Y M(-P) • f(G/ P) = Ya(G) = 
Pen(G) 

for all G with |F(G)| > n = f(K\), where the first equality holds by ( |2.8| ) and the second because for 
G with |F(G)| > n there are no injective functions a : V(G) — > [n\. 

It is quite surprising that these trivial conditions are sufficient to guarantee that / is a partition 
function. To see that they are, let / : Q — > C be a multiplicative graph invariant such that ( |2.10[ ) holds 
for all G G G with \V{G)\ > f(Ki). 

We first show that n := f{K{) is a non-negative integer. Let k := [|,f (iTi)|] , and let 1^ be the graph 
with vertex set [k] and no edges. (Hence = if k = 0.) Suppose that n is not a non-negative integer. 
Then k > \n\, and by ( |2.10[ ), the multiplicativity of /, and ( |2.9[ ) we have 



0= Yl KP)-f(h/P)= Yl v(P)-n lPl =n-(n-l)---(n-k + l) / 0. 

p e n(/ fc ) Pen([fc]) 

This is a contradiction. 

A quantum graph is a formal linear combination of graphs with coefficients from C, that is, an 
expression Yli=i a iGi, where £ > and a-i G C and G-i G Q for all i G [£]. The class of all quantum 
graphs is denoted by QQ. The quantum graphs obviously form a vector space over C, and by extending 
the product "disjoint union" linearly from Q to QQ, we turn this vector space into an algebra. We also 
extend the function / linearly from Q to QQ. Observe that / is an algebra homomorphism from QQ to 
C, because it is multiplicative. 

For all i, j G [n] we let X^jy be a variable, and we let X be the tuple of all these variables ordered 
lexicographically. Furthermore, we let X be the n x n matrix with Xij := Xjj := X^jy. We view X 



2.2 Which functions are partition functions ? 
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as a matrix over the polynomial ring C[X]. We extend the function Zx '■ G — > C[X] linearly from Q to 
QQ. Then Zx is an algebra homomorphism. It is not too hard to show that the image Zx{QG) consists 
precisely of all polynomials in C[X] that are invariant under all permutations Xuj\ h-» Xt^U),^)} °f 



the variables for permutations it of [re]. Using (2.8 1 and other properties of the Mobius inversion, it can 



be shown that the kernel of Zx is contained in the kernel of /. This implies that there is an algebra 
homomorphism / from the image Zx{QQ) to C such that / = / o Zx- 
Then 

/:= {peZx(QG)\M = 0} 

is an ideal in the subalgebra Z X (QQ) Q C[X]. We claim that the polynomials in / have a common zero. 
Suppose not. Let I' be the ideal generated by I in C[X]. Then by Hilbert's Nullstellensatz it holds that 
1 G I'. Using the fact that the the subalgebra Zx(QQ) consists precisely of all polynomials in C[X] 
that are invariant under all permutations Xuj\ i-> X{n(i),*(j)} f° r n £ S n , one can show that actually 
lei But then 

Q = f(l) = f(Z x (K )) = f(K ) = l, 

which is a contradiction. 

Thus the polynomials in / have a common zero. Let A = (Ar»j) | i,j £ [re]) be such a zero, 
and let A G C nxn be the corresponding symmetric matrix. Observe that for each graph G it holds that 
Z X (G) - f(G) € I, because f(Z x (G)) - /(/(G)) = /(G) - /(G) = 0. Hence Z A (G) - /(G) = 
and thus f(G) = Za(G). This completes our proof sketch. □ 



Even though Theorems |2. 1 1 1 and 2.12 look quite different, it is not hard to derive Theorem 2.11 from 



Theorem |2.12[ In the remainder of this subsection, we sketch how this is done. 



Proof of Theorem 2.11 (sketch). We have already proved the forward direction. For the backward 



direction, let / : Q — > R be a multiplicative graph invariant such that all connection matrices of / are 



positive semidefinite. We first prove that / satisfies condition ( |2.10| ): Let n := f(K\) (we do not know 
yet that re is an integer, but it will turn out to be), and let k > re be a non-negative integer. Let 1^ be the 
graph with vertex set [k] and no edges, and for every partition P of [k], define <pp : [k] — > V{Ik/P) to 
be the canonical projection, that is, 4>p(i) is the class of i in the partition P. Then (Ik/P, <ftp) £ Qk- As 
the fcth connection matrix M(f, k) is positive semidefinite, we have 

Y, KP) ■ t*(Q) ■ f{(h/P 4>p) ■ (h/Q, ^q)) > o. 

p,Qen([fc]) 

For partitions P,QG II([A;]), let P V Q be the least upper bound of P and Q in the partially ordered set 
II([A;]), and note that (Ik/ P,4>p) ■ {Ik/Q,<t>Q) = (Ik/ P V Q, 4>p\jq)- Hence by the multiplicativity of 
/ we have f((Ik/P, 4>p) ■ (Ik/Q, </>q)) = rel Pv< ^L A calculation similar to the one that leads to ( 2.9 1 



shows that £p.Qen[fc] n(P) ■ «(Q) ■ Xl PvQ l = X ■ (X - 1) • • • (X - k + 1). Hence for all non-negative 
integers k > re we have 

0< Yl v(P)-riQ)-f[(Ik/P,<t>p)-(Ik/Q,<i>Q))=n-(n-l)---(n-k + l). (2.11) 
p,Qen([fc]) 

This is only possible if re is a non-negative integer, in which case equality holds. 

Now let G be an arbitrary graph with k := | V(G)| > re. Without loss of generality we may assume 
that V(G) = [k). Let ip be the identity on [k]. Consider the (|n([A;])| + 1) x (|n([A;])| + 1) -principal 
submatrix Mq of M(f, k) with rows and columns indexed by the /c-labeled graphs (Ik/P, 4>p) f° r ai l 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



P G n([/c]) and (G,ip). For every x G M., let v x be the column vector with entries fi(P) for all 
P G n([A;]) followed by x. By the positive semi-definiteness of M(f, k), we have 

< vjMoVa, 

= KP)-KQ)-f{(h/P,4>p)-(h/Q,4> Q )) 
p,Qen([fc]) 

+ 2x- J] M^)-/((4/ J P,0p)-(G,^))+x 2 -/((G,V)-(G,^)) 
p e n([fc]) 

= 2x- Yl KP)-f((h/P,<pp)-(G,iP))+x 2 -f{(G,iP)-(G,^) 
pen([fc]) 



This is only possible for all x G K if ^Pen(M) M-P) ' f{(h/P, <j>p) • (G, tp)) = 0. Note that for every 
P G U([k}) it holds that {I k /P,<j> P ) ■ (G,%/>) = (G/P,<f> P ). Thus 

Y KP)'f(G/P) = 0. 
Pen(G) 



If follows from Theorem 2.12 that / = Za for some matrix A G C nxn . 

For every i > 0, let Fi be the graph with two vertices and £ parallel edges between these vertices. 
Exploiting the positive semi definiteness of the principal submatrix of M(f, 2) with rows and columns 
indexed by the graphs Fg for < £ < n(n + l)/2, it is not hard to show that the matrix A is real. □ 



2.3 Generalizations 
Vertex weights 

We have mentioned that a partition function Za can be viewed as mapping a graph G to a weighted sum 
of homomorphisms to the edge-weighted graph represented by the matrix A. We may also put weights 
on the vertices of a graph. It will be most convenient to represent vertex weights on a graph H by a 
diagonal matrix D = (D V)W ) v we yr H \, where D v>v is the weight of vertex v and D V)W = if v ^ w. 
Then we define the weight of a mapping a from G to H to be the product of the weights of the images 
of the edges and the images of the vertices. More abstractly, for every symmetric matrix A G § nxn and 
diagonal matrix D G § nxn we define a function Za,d '■ Q — > § by 

Za,d(G) := Y, II A cr{v),a{w) ' \\ D a(v),a(v)- 

cr:V(G)-¥[k] vweE(G) v£V(G) 

We call Za,d a partition functions with vertex weights over S. ([Freedman, Lovasz, and Schrij ver MFLS07H 
use the term homomorphism functions.) The vertex weights enable us to get rid of the constant factors 
in some of the earlier examples and give smoother formulations. For example: 



Example 2.13. Recall that by Example 2.9 the number F(G, k) of nowhere-zero k-flows of an N -vertex 
graph G is k~ N -ZA{G)forthe kxk matrix A with diagonal entries (k — 1) and off -diagonal entries —1. 
Let D be the (k x k)-diagonal matrix with entries Da := 1 /kfor all i G [k]. Then F(G, k) = Za,d(G). 



Freedman, Lovasz, and Schrijver gave an algebraic characterization of the class of partition functions 
with non-negative vertex weights over the reals, again using connection matrices. However, they only 



2.3 Generalizations 
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consider graph invariants defined on the class Q' of graphs without loops (but with parallel edges). For 
a graph invariant / : Q' — > M, we define the kth connection matrix M(f, k) as for graph invariants 
defined on Q, except that we omit all rows and columns indexed by graphs with loops. We call a matrix 
A G M nxn non-negative if all its entries are non-negative. 

Theorem 2.14 (Freedman, Lovasz, and Schrijver |FLS07|). Let f : Q' — > R be a graph invariant. 
Then the following two statements are equivalent: 

1. There are a symmetric matrix A G M. nxn and a non-negative diagonal matrix D G M nxn such 
that f(G) = Z AyD (G)forall G G Q' . 

2. There is a q > such that for all k > the matrix M(/, k) is positive semidefinite and has row 
rank at most q k . 



Despite the obvious similarity of this theorem with Theorem |2.1 1| the known proofs of the two theorems 
are quite different. 



Asymmetric matrices and directed graphs 

Of course we can also count homomorphism between directed graphs. We denote the class of all directed 
graphs by V; as undirected graphs we allow directed graphs to have loops and parallel edges. For every 
directed graph H, we define the homomorphism counting function Zh '■ T> — > Z by letting Zjj(D) be 
the number of homomorphisms from D to H. For every square matrix A G S nx ™ we define the function 

Z A : V -> S by 

Za{D) := II A <r(v),cr(w)- 

a:V(D)-t[n] (v,w)&E(D) 

Note that if A is a symmetric matrix, then Z A {D) = Z A (G), where G is the underlying undirected 
graph of G. 



Hypergraphs and relational structures 

Recall that a hypergraph is a pair H = (V, E) where V is a finite set and E Q 2 V . Elements of V 
are called vertices, elements of E hyperedges. A hypergraph is r-uniform, for some r > 1, if all its 
hyperedges have cardinality r. Thus a 2-uniform hypergraph is just a simple graph. A homomorphism 
from a hypergraph G to a hypergraph if is a mapping h : V(G) — > V(H) such that h(e) G E(H) for all 
e G E(H). Of course, for e = {vi, . . . , v r } we let h(e) = {h(v\), . . . , fr(tv)}Q Hypergraph homomor- 
phism counting functions and partition functions have only been considered for uniform hypergraphs. 
The class of all r-uniform hypergraphs is denoted by fi r . For all H G T-L r , we define the homomorphism 
counting function Zh ■ rl r — > Z in the obvious way. The natural generalization of partition functions 
to r-uniform hypergraphs is defined by symmetric functions A : [n} r — > S. We define Z A by 

Za(G):= Y, II A ^ e ^ 
a:V(G)->[n] e£E(G) 

where for e = {v\, . . . , v r } we let A(a(e)) := A(a(vi), . . . , o~(v r )). This is well-defined because A is 
symmetric. 



2 There is no need to define hypergraph homomorphisms as mappings from V(G) U E{G) to V{H) U E(H) because we 
do not allow parallel hyperedges. 
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2 FROM HOMOMORPHISMS TO PARTITION FUNCTIONS 



Let us finally consider homomorphisms between relational structures. A (relational) vocabulary a 
is a set of relation symbols; each relation symbol comes with a prescribed finite arity. A a -structure B 
consist of a set V(B), which we call the universe or vertex set of B, and for each r-ary relation symbol 
S 6 a an r-ary relation 5(B) C V(B) r . Here we assume all structures to be finite, that is, to have 
a finite universe and vocabulary. For example, a simple graph may be viewed as an {E} -structure G, 
where E is a binary relation symbol and E(G) is irreflexive and symmetric. An r-uniform hypergraph 
H may be viewed as an {E r } -structure, where E r is an r-ary relation symbol and E r (H) is symmetric 
and contains only tuples of pairwise distinct elements. For each vocabulary a, we denote the class 
of all ex-structures by S a . For each cr-structure B we define the homomorphism counting functions 
Zb '■ S a — > Z in the usual way. 

To generalize partition functions, we consider weighted structures. For a ring § and a vocabulary 
er, an ^-weighted a-structure B consists of a finite set V(B) and for each r-ary R £ a a mapping 
R(S) : V{B) r -> §. Then we define Z B : S a -»• § by 

^(A):= EII II R(S)(a( ai ),...,a(a r )). 

a:V(G)-+{n] Recr (ai,...,a r )eV(A) r , 
where r is the arity of i? 

Note that this does not only generalize plain partition functions on graphs, but also partition functions 
with vertex weights, because we may view a graph with weights on the vertices and edges as a weighted 
{E, P}-structure, where E is a binary and P a unary relation symbol. 

It is a well-known observation due to Feder and Vardi [FV98] that constraint satisfaction problems 
may be viewed as homomorphism problems between relational structures and vice versa. Thus counting 
homomorphisms between relational structures correspond to counting solutions to constraint satisfaction 
problems. 



Edge Models 

Partition functions and all their generalizations considered so far are weighted sums over mappings 
defined on the vertex set of the graph or structure. In the context of statistical physics, they are sometimes 
called vertex models. There is also a notion of edge model. An edge model is given by a function 
F : N n — > 8. Let G = (V, E) be a graph. Given a "configuration" r : E — >• [n], for every vertex v we 
let t(r, v) := (ti,..., t n ), where ti is the number of edges e incident with v such that r(e) = i. We 
define a function Zp : Q — > § by 

Zf{G):= E \{F(t^v)). 

r:E-¥[n] f£V 

Szegedy [Sze07] gave a characterization of the graph invariants over the reals expressible by edge mod- 



els that is similar to the characterizations of partition functions (vertex models) given in Theorems 2.11 
andl2~T4l 

Example 2.15. Let F : N 2 — > Z be defined by F(i,j) := 1 if j = 1 and F(i,j) = otherwise. 
Let G = (V, E) be a graph. Observe that for every r : E(G) — > {1,2} and every v £ V it holds 
that F(t(T,v)) = 1 if and only if there is exactly one edge in r _1 (2) that is incident with v. Hence 
Ylv&v ^ ? (*( r i v )) = ^ if t^ 1 (2) is a perfect matching of G and W veV F{t(T,v)^ = otherwise. It 
follows that Zp(G) is the number of perfect matchings of G. 

It can be proved that the function f = Zp counting perfect matchings is not a partition function, 
because the connection matrix M(f, 1) is not positive semi-definite IF LS07\I . 
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Related to the edge models is a class of functions based on Valiant's holographic algorithms [Val08 |: 
so-called holant functions which have been introduced by Cai, Lu, and Xia [CLX09]. A holant function 
is given by a signature grid Q = (G, J 7 , it) where G = (V, E) is a graph and, for some n 6 N, the set 
T contains functions / : [n] a f — > §, each of some arity a,f. Further it maps vertices v G V to functions 
f v G T such that d(v) = af v . Let, for some vertex v, denote E(v) the set of edges incident to v. The 
holant function over this signature grid is then defined as 

Holant^ = II fv{T\ E (v))- 

r:E-¥[n] v£V 

Note that we assume here implicitly that E has some ordering, therefore f v {T~\E(v)) is well-defined. 
Edge models are a special case of this framework, since Zp{G) = Holant^ for the signature grid which 
satisfies f v {T\E(v)) = F{t{ T i v )) f° r ai l v E V. More generally it can be shown that partition functions 
on relational structures can be captured by holant functions. 



3 Complexity 

Partition functions tend to be hard to compute. More precisely, they tend to be hard for the complexity 
class #P introduced by Valiant in [Val79], which may be viewed as the "counting analogue" of NP. A 
counting problem C (that is, a function with values in the non-negative integers) belongs to #P if and 
only if there is a nondeterministic polynomial time algorithm A such that for every instance x of C it 
holds that C(x) is the number of accepting computation paths of A on input x. There is a theory of 
reducibility and #P-completeness much like the theory of NP-completeness. (The preferred reductions 
in the complexity theory of counting problems are Turing reductions. We exclusively work with Turing 
reductions in this article .) Most NP-complete decision problems have natural #P-complete counting 
problems associated with them. For example, the problem of counting the number of independent sets of 
a graph and the problem of counting the number of 3-colorings of a graph are both #P-complete. There 
is a counting analogue of a well-known Theorem due to Ladner [Lad75] stating that there are counting 
problems in #P that are neither #P-complete nor in FP (the class of all counting, or more generally 
functional problems solvable in polynomial time); indeed the counting complexity classes between FP 
and #P form a dense partial order. 

Even though the class of partition functions is very rich, it turns out that partition functions and also 



their various generalizations discussed in Section 2.3 exhibit a complexity theoretic dichotomy: Some 



partition functions can be computed in polynomial time, most are #P-hard, but there are no partition 
functions of intermediate complexity. A first dichotomy theorem for homomorphism counting functions 
of graphs was obtained by Dyer and Greenhill: 

Theorem 3.1 (Dyer and Greenhill |DG00|). Let H be a graph without parallel edges. Then Zh is 
computable in polynomial time if each connected component of H is either a complete graph with a 
loop at every vertex or a complete bipartite graph^\ Otherwise, Zh is ^-complete. 

When thinking about the complexity of partition functions over the real or complex numbers, we face 
the problem of which model of computation to use. There are different models which lead to different 
complexity classes. To avoid such issues, we restrict our attention to algebraic numbers, which can 



3 We count the empty graph and the graph K\ with one vertex and no edges as bipartite graphs. 
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3 COMPLEXITY 



be represented in the standard bit modelj^] We will discuss the representation of and computation with 



algebraic numbers in Section 7.3 The fields of real and complex algebraic numbers are denoted by Ra 
and Ca, respectively. In general, partition functions are no counting functions (in the sense that their 
values are no integers) and hence they do not belong to the complexity class #P. For that reason, in the 
following results we only state #P-hardness and not completeness. As an easy upper bound, we note 
that partition functions over Ra and Ca belong to the complexity class FP* P of all functional problems 
that can be solved by a polynomial time algorithm with an oracle to a problem in #P. 

Let us turn to partition functions over the reals. We first observe that if A £ R^ Xn 1S a symmetric 
matrix of row rank 1, then Za is easily computable in polynomial time. Indeed, write A = a T a for a 
(row) vector a = (a±, . . . , a n ) G R^. Then for every graph G = (V, E), 

n 

Za{G)= £ Ua a{v) a aiw) = £ n^f^IIE^- (3-D 

a:V^[n] vw&E a:V-+{n] v£V vGV i=l 

Here deg(v) denotes the degree of a vertex v. The last term in Q3.1) , which only involves a polynomial 
number of arithmetic operations, can easily be evaluated in polynomial time. Thus partition functions 
of rank-1 matrices are easy to compute. It turns out that all easy partition functions of non-negative real 
matrices are based on a "rank-1" condition. 

Let us call a matrix A bipartite if its underlying graph G(A) is bipartite. After suitable permuting 
rows and columns, a bipartite matrix has the form 

B 
B T 

where we call B and B T the blocks of A. A similar argument as the one above shows that if B has 
row rank 1 then Za is computable in polynomial time. Note that we only need to compute Za{G) for 
bipartite G, because Za(G) = for non-bipartite G. 

The connected components of a matrix A are the principal submatrices corresponding to the con- 
nected components of the underlying graph G{A). Note that if A is a matrix with connected components 
A\, . . . , A m , then for every connected graph G it holds that Za{G) = Y^=x ^Aj{G). By the multi- 
plicativity of partition functions, for a graph G with connected components G\, . . . , Gg, we thus have 
Za{G) = Y[i=i Sj=i (Gi). This reduces the computation of Za to the computation of the Za 3 - 

Theorem 3.2 (Bulatov and Grohe |BG05 |). Let A 6 R2 X ™ b e a symmetric non-negative matrix. Then 
Za is computable in polynomial time if all components of A are either of row rank 1 or bipartite and 
with blocks of row rank 1. Otherwise, Za is #P-hard. 



Note that the theorem is consistent with Theorem 3.1 the special case for 0-1 -matrices. We have already 



proved that Za is computable in polynomial time if all components of A are either of row rank 1 or 
bipartite and with blocks of row rank 1. The much harder proof that all other cases are #P-hard will be 
given in Sections [4]-[9] of this paper. As a by-product, we also obtain a proof of Theorem |3.1| 

The rest of this section is a survey of further dichotomy results. We will not prove them. The 
following two examples show that if we admit negative entries in our matrices, then the rank-1 condition 
is no longer sufficient to explain tractability. 



4 In |BG05 1, the complexity classification of partition functions for non-negative real matrices (Theorem 3.2 of this article) 



was stated for arbitrary real matrices in one of the standard models of real number computation. However, the proof is faulty 



and can only be made to work for real algebraic numbers (cf. Section 7.3 of this article) 
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Example 3.3. Consider the matrix B 



first introduced in Example 



1 1 

1 -1 

this matrix is 2, yet we shall prove that Zb is computable in polynomial time. Let 



2.6 

G 



The row rank of 
(V, E). Then 



Z B (G) 



E 

o-:V-+[2] vw&E 



n 



a:V-^[2] vw&E 
= ^2 (-l)^^"eB< J W- CT ( w ). 
<r:V-¥[2] 

Hence Zb(G) is 2 N minus twice the number of mappings a : V — > {0, 1} such that YlvweE a ( v ) ' o~(w) 
is odd. Here N is the number of vertices of G. Thus to compute Zb, we need to determine the number 
of solutions of the quadratic equation 

vw£E 

over the 2-element field ¥2. The number of solutions of a quadratic equation over ¥2 or any other finite 
field can be computed in polynomial time. This follows easily from standard normal forms for quadratic 
equations (see, for example, [LN97], Section 6.2). 

" xn and B G S kxl is the m ■ k x n ■ £ matrix 
A ltn -B 



A ■ R 



Example 3.4. The tensor product of two matrices A G §' 

(A ltX B 
A®B:= : 

\A m ,i ■ B 

It is easy to see that for all square matrices A, B and all graphs G it holds that 

Za®b{G) = Za{G) ■ Z B {G). 

We can thus use the tensor product to construct new matrices with polynomial time computable partition 
functions. For example, the following three 4 x A-matrices have polynomial time computable partition 
functions: 



( 


1 


1 


1 


1 


\ 


/ 


-1 


-1 


1 


1 \ 




( 1 


2 


1 


2 \ 




1 


-1 


1 


-1 






-1 


1 


1 


-1 




2 


4 


2 


4 




1 


1 


-1 


-1 






1 


1 


-1 


-1 




1 


2 


-1 


-2 


V 


1 


-1 


-1 


1 


/ 


V 


1 


-1 


-1 


1 J 




\2 


4 


-2 


-4 ) 



Very roughly, all symmetric real matrices with a polynomial time computable partition function 
can be formed from matrices of rank 1 and matrices associated with quadratic equations over F2 (in a 
way similar to the matrix B of Example 3.3 ) by tensor products and similar constructions. The precise 
characterization of such matrices is very complicated, and it makes little sense to state it explicitly. In 
the following, we say that a class T of functions exhibits an FP - j^V -dichotomy if all functions in F 
are either in FP or #P-hard. We say that F exhibits an effective FP - #P-dichotomy if in addition it is 
decidable if a given function in F, represented for example by a matrix, is in FP or #P-hard. 

Theorem 3.5 (Goldberg, Grohe, Jerrum, Thurley HGGJTlOll ). The class of partition functions of sym- 
metric matrices over the reals exhibits an effective FP - #P- dichotomy. 



There are two natural ways to generalize this result to complex matrices. Symmetric complex matrices 
where studied by Cai, Chen, Lu HCCLlOal . and Hermitian matrices by Thurley [Thu09]. 
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4 AN ITINERARY OF THE PROOF OF THEOREM ?? 



Theorem 3.6 (Cai, Chen, Lu BCCLlOal . Thurley |Thu09|). 1. The class of partition functions of 
symmetric matrices over the complex numbers exhibits an effective FP - #P- dichotomy. 

2. The class of partition functions of Hermitian matrices exhibits an effective FP - #P- dichotomy. 
Beyond Hermitian matrices, partition functions of arbitrary, not necessarily symmetric matrices are 



much more difficult to handle. An FP - #P-dichotomy follows from Bulatov's Theorem 3.8 below, but 
it is difficult to understand how this dichotomy classifies directed graphs. Dyer, Goldberg and Paterson 
[ DGP07 ] proved an effective FP - #P-dichotomy for the class of homomorphism counting functions Zh 
for directed acyclic graphs H; it is based on a complicated "rank-1" condition. A very recent result by 
Cai and Chen HCC101 establishes an effective dichotomy for partition functions Za on non-negative real- 
valued (i.e. not necessarily symmetric) matrices A. The complexity of hypergraph partition functions 
was studied by Dyer, Goldberg, and Jerrum [DGJ08 ], who proved the following theorem: 

Theorem 3.7 (Dyer, Goldberg, and Jerrum |DGJ08|). For every r, the class of functions Z a from the 
class 7~L r of r -uniform hypergraphs defined by non-negative symmetric functions A : [n] r — > Ma exhibits 
an effective FP - jfi' -dichotomy. 

Let us finally turn to homomorphism counting functions for arbitrary relational structures, or equiva- 
lently solution counting functions for constraint satisfaction problems. Creignou and Hermann [CH96] 
proved a dichotomy for the Boolean case, that is, for homomorphism counting functions of relational 
structures with just two elements. Dichotomies for the weighted Boolean case were proved in [DGJ09] 
for non-negative real weights, in ||BDG + 09l for arbitrary real weights, and in [CLX09] for complex 
weights. Briquel and Koiran [BK09] study the problem in an algebraic computation model. The general 
(unweighted) case was settled by Bulatov: 

Theorem 3.8 (Bulatov |Bul08|). The class of homomorphism counting functions for relational struc- 
tures exhibits an FP - #P- dichotomy. 

Bulatov's proof uses deep results from universal algebra, and for some time it was unclear whether his 
dichotomy is effective. Dyer and Richerby [DRIOb] gave an alternative proof which avoids much of 
the universal algebra machinery, and they could show that the dichotomy is decidable in NP (DRar, 
IDRlOaL Extensions to the weighted case for nonnegative weights in Q have been given by Bulatov et 
al. ||BDG + I0t . and for nonnegative algebraic weights by Cai et al. BCCLlObB . 

In the context of holant functions complexity dichotomies have been obtained by Cai, Lu, and Xia 
BCLX081 ICLX09 1 . Besides the result stated above, they considered holant functions on signature grids 
£1 = (G, F, it), given that J 7 is a set of symmetric functions satisfying certain additional conditions. 
A first a dichotomy [CLX08] is given for any of boolean symmetric functions F on planar bipartite 
2, 3-regular graphs. In [CLX09] dichotomies are presented, assuming that the class F contains certain 
unary functions. 



4 An Itinerary of the Proof of Theorem 3.2 



From now on, we will work exclusively with partition functions defined on matrices with entries in one 
of the rings QPH > Q> 7l\X\ , Z, and we always let § denote one of these rings. For technical reasons, 
we will always assume that numbers in Ma are given in standard representation in some algebraic 
extension field Q(9). That is, we consider numbers in Q(8) as vectors in a (i-dimensional Q-vectorspace, 
where d is the degree of Q(9) over Q. It is well-known that for any finite set of numbers from Ma we 



17 



can compute a 6 which constitutes the corresponding extension field (cf. |Coh93] p. 181). For further 
details see also the treatment of this issue in [DGJ08i rThu09i 

By deg(/) we denote the degree of a polynomial /. For two problems P and Q we use P < Q to 
denote that P is polynomial time Turing reducible to Q. Further, we write P = Q to denote that P < Q 
and Q < P holds. 

An m x n matrix A is decomposable, if there are non-empty index sets / C [m], J C [n] with 
(I, J) [m] x [n] such that Aij = for all (i, j) & I X J and all (i, j) & I X J, where J := [m] \ I 
and J := [ra] \ J. A matrix is indecomposable if it is not decomposable. A block of A is a maximal 
indecomposable submatrix. Let A be an m x m matrix and G := its underlying graph. Note that 
every connected component of G that is not bipartite corresponds to a block of A, and every connected 

B 
B T 

The proof of Theorem |3 . 2| falls into two parts, corresponding to the following two lemmas. 

Lemma 4.1 (Polynomial Time Solvable Cases). Let A G § mxm b e symmetric. If each block of A has 
row rank at most 1 then EVAL(A) is polynomial time computable. 

We leave the proof of this lemma as an exercise for the reader. The essential ideas of its proof were 
explained in Section [3] before the statement of Theorem 3.2 



component that is bipartite corresponds to two blocks B, B T arranged as in 



Lemma 4.2 (#P-hard Cases). Let A E § mxm \j e symmetric and non-negative. If A contains a block of 
row rank at least 2 then EVAL(j4) is #V-hard. 



Theorem 3.2 directly follows from Lemmas 4.1 and 4.2 The remainder of this paper is devoted to a 



proof of Lemma 4.2 



For technical reasons, we need to introduce an extended version of partition functions with vertex 
weights. Several different restricted flavors of these will be used in the proof to come. Let A € M™ xm 
be a symmetric matrix and D G M™ xm a diagonal matrix. Let G = (V, E) be some given graph. Recall 
that a configuration is a mapping a : V — > [m] which assigns a spin to every vertex of G. By contrast, 
a pinning of (vertices of) G with respect to A is a mapping <p : W — > [m] for some subset W C V. 
Whenever the context is clear, we simply speak of pinnings and configurations without mentioning the 
matrices A,D and the graph G explicitly. We define the partition function on G and <ft by 

(f>Ca:V-^[rn] uveE ■ueV\dom(0) 

where dom(</>) denotes the domain of 4>. Note that the sum is over all configurations a : V — > [m] 
which extend the fixed given pinning (p. In the presence of a pinning <ft we denote the weight of the 
configuration a by the following term 

(v),a(v) • 

uvGE ueV\dom(<£) 

Note that, for technical reasons, the terms D a / V \ a / V \ for v G dom((/>) are excluded from this weight. 
Whenever cj> is empty in the sense that dom(</>) = then we say that it is trivial. In this case, its 
appearance in the above expression is vacuous. This is analogously true for D if it is the identity 
matrix. In either of these cases we omit the terms D (0, respectively) in the expression. For example if 
dom((/>) = and D = I m , then 

z a(G) = ^2 IT A v(u),a(v)- 
cr:V—t\m] uv£E 
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4 AN ITINERARY OF THE PROOF OF THEOREM ?? 



The definitions of Za,d(G) and Za(<(>, G) are analogous. We define EVAL pm (j4, D) as the computa- 
tional problem of computing Za,d(4>, G) on input 4>, G. Similarly EVAL(j4, D) restricts the inputs to 
empty pinnings and EVAL pin (^4) denotes the problem where D is the identity matrix. 



We will now explain the overall structure of the proof of Lemma 4.2 In a first step we will see how 



we can augment our capabilities so as to fix some vertices of the input graphs, without changing the 



complexity of the problems under consideration (cf. Lemma 4.3 ). Then in the General Conditioning 



Lemma 4.4 we will show, that we can reduce the abundance of non-negative matrices to certain well- 



structured cases. From these we will show in two steps (Lemmas 4.5 and 4.6 ) how to obtain #P-hardness. 

Lemma 4.3 (Pinning Lemma). Let A G IR™*" 1 be a symmetric non-negative matrix. Then 

EVAL p/n (^) = EVAL(A). 

A 1-cell in a matrix A G § mxm i s a submatrix Ajj such that Aij = 1 for all (i, j) G I x J and A+j ^ 1 
for all G (I x J) U (I x J). For a number or an indeterminate X an X -matrix is a matrix whose 
entries are powers of X. 

General Conditions. For a matrix A G § mxm we define conditions 

(A) A is symmetric positive and has rankvl > 2. 

(B) A is an X-matrix for an indeterminate X. 

(C) There is a k > 2, numbers 1 = uiq < . . . < m& = m + 1 and Ii = [raj_i, mi — 1] for all i G [fc] 

such that Aj i j i is a 1-cell for every i G [fc — 1]. The matrix Af fc / fc may or may not be a 1-cell. 
Furthermore, all 1-entries are contained in one of these 1-cells. 

Lemma 4.4 (General Conditioning Lemma). Let A G § mxm be a non-negative symmetric matrix 
which contains a block of rank at least 2. Then there is a 7L\X\-matrix A satisfying conditions (A)-(C) 
such that 

EVAL p/ V) < EVAL p/n (,4). 

If a matrix A satisfies the General Conditions two different cases can occur which we will treat separately 
in the following. The first case is the existence of at least two 1-cells. 

Lemma 4.5 (Two 1-Cell Lemma). Let A G Z[X] mxm be a positive symmetric matrix containing at 
least two 1-cells. Then EVAL p/ "(,4) is #P-hard. 

The second case is then the existence of only one 1-cell. The proof of this case is much more involved 
than the first one. 

Lemma 4.6 (Single 1-Cell Lemma). Let A G Z[X] mxm be a matrix which satisfies conditions (A) - 
(C) and has exactly one 1-cell. Then EVAL f " n (yl) is j^P-hard. 

Once these results have been derived, it will be easy to prove the main result. 



Proof (of Lemma 4.2). Let A G MY 1 m be a non-negative symmetric matrix which contains a block of 



rank at least 2. By the General Conditioning Lemma 4.4 there is a matrix A satisfying conditions (A) 



(C) such that EVALP m (^') < EVAL pm (A). If A contains a single 1-cell then EVALP m (A') is #P-hard 



by Lemma |4.6| Otherwise, it is #P-hard by Lemma |4.5| In both cases this proves #P-hardness of 
EVAL(A) by means of Lemma 4.3 n 
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5 Evaluation and counting 

Let Abe an m x m matrix. Then, for every q G Q we define the matrix by 



A {i) = { ( A ij) q , XAij^O 
tJ \ , otherwise. 



The following lemma provides two basic reductions which form basic building blocks of many hardness 
proofs. 

Lemma 5.1. Let A G § mxm the following is true for every p G N 
(p-thickening) EVAL pin (A^) < EVAL pin (A). 
(p-stretching) EVAL pin (Ap) < EVAL pin (A). 

Proof. Let G = (V, E) be a graph and <f> a pinning. Let the p-thickening of G be the graph obtained 
from G by replacing each edge by p many parallel edges. The p-stretching G p of G is obtained from G 
by replacing each edge by a path of length p. The reductions follow from the (easily verifiable) identities 

Z A(p) (<f>, G) = Z A (4>, G(p)) and Z AP (<P, G) = Z A (<f>, G p ). " □ 

Lemma 5.2. Let A G § mx ™ be non-negative. If A contains a block of rank > 2 then so does AA T . 

Proof. Let Ajj be a block of rank at least 2 in A and define A' = AA T . We claim that A' n contains a 
block of rank at least 2, which clearly implies the statement of the lemma. 

For G / let Ai^ and A^i * be two linearly independent rows. Since Ajj is a block, there are 
indices i = ii, . . . ,ik = i' in / and ji,.. -jk-i G J such that Ai v j v / and Ai v+1 j v / for all 
v G [fe — 1]. Linear independence of Aj,^ and implies that there is a ^ G [k — 1] such that also 
^4j Mi * and are linearly independent. We claim that the submatrix 

is a witness for the existence of a block of rank at least two in A'. To see this note first that, by definition, 
all entries of this submatrix are positive, it thus remains to show that this submatrix has rank 2. Assume, 
for contradiction, that it has zero determinant, that is 

(^4« M ,*, ^4j M ,*)(^4j M+ i,*, ^4j M+ i,*) = (^ip+i,*!^,*) • 

The Cauchy Schwarz inequality therefore implies linear dependence of ^4^,* and ^4j M+li *. Contradic- 
tion, n 

For A an m x m matrix and ir : [m] — > [to] a permutation, define A wn by 

(A-Tr)ij = A n(i)n(j) for all i, j G [m]. 

The Permutability Principle states that for any evaluation problem on some matrix A, we may assume 
any simultaneous permutation of the rows and columns of this matrix. Its proof is straightforward. 

Lemma 5.3 (Permutability Principle). Let A,D G § mxm and it : [m] — > [m] a permutation. Then 
E\AE pin {A,D)=E\AE pin {A^,D^). 
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5 EVALUATION AND COUNTING 



We will make extensive use of interpolation; the following simple lemma is one instance. 

Lemma 5.4. For some fixed 6 E Ma let x\, . . . , x n E Q(#) be pairwise different and non-negative 
reals. Let b±, . . . ,b n E Q(0) &e arbitrary such that 

n 

bj = °i x i f or Q H 3 £ [ n \ ■ 
i=l 

r/jen the coefficients c\ , . . . , c n are uniquely determined and can be computed in polynomial time. 
5.1 The Equivalence of EVAL pin (A) and COUNT pin (A) 

Let A E S mxm be a matrix, G = (V, E) a graph and <fi a pinning. We define a set of potential weights 

Wa(G) := < JJ ^ | ^ 771^ = 1^1, and mij > 0, foralli, j E [m] > . (5.1) 
\i,je[m] i,ie[m] I 



For every w£§ define the value 



N A (G,<t>,w) :-- 



a : V ->■ [m] 



By COUNT pin (A) we denote the problem of computing Na{G, 0, w) for given G, <p and to E S. In 
analogy to the evaluation problems we write COUNT(A) for the subproblem restricted to instances 
with trivial pinnings. It turns out that these problems are computationally equivalent to the evaluation 
problems of partition functions. 

Lemma 5.5. For every matrix A E § mxm we have 

EVAL pin (A) = COUNT p/ "(A) and EVAL(A) = COUNT(A). 
Proof. Let G = (V, E) be a graph and (f> a pinning. We have 

Z A ((j),G) = ^ II A criu)a{v) = ^2 w ■ N A(G,(f),w). 
<t>Co-.V->[m]uv£E w€_W A (G) 

As the cardinality of W A (G) is polynomial in the size of G this proves the reducibilities 

EVAL pin (A) < COUNT pin (A) and EVAL(A) < COUNT(A). 

For the backward direction let G^ denote the graph obtained from G by replacing each edge by t 
parallel edges. We have 

Z A {^G^)= fil^(-M-)) = E w l -N A {GAM- 

<j>Ca:V-^[m] \uv£E J weW A (G) 

Using an EVAL pin (A) oracle, we can evaluate this for t = 1, . . . , |WU(G)|- Therefore, if S is one of 



x , Q, Z, the values N A (G, (f>, w) can be recovered in polynomial time by Lemma 5.4 



5.2 Dealing with Vertex Weights 
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If § is one of Z[X] or Q[X], then let f t (X) denote Z A ((f>, G (t) ). We obtain the following system of 
equations, for t = 1, . . . , \ Wa 

f t (X)= w\X)-N A {G,4>,w). (5.2) 

w£\V A (G) 

Let det(X) denote the determinant of this system of equations. W.l.o.g., all the w(X) G W A (G) are 
non-zero polynomials, det(X), as it is a Vandermonde determinant in these w(X), is itself a non-zero 
polynomial of the form 

det(X)= w(X). ft (w(X)-w'(X)). 

Let<5 = 1 + max{degu>(X) | w G W^G)} and observe that each w G VVa(G) has at most <5 — 1 roots. 
Further, each of the terms w(X) — w'{X) of det(X) has degree at most 5 — 1 and thus the degree of 
det(X) is strictly smaller than 5' := (\ Wa ^\ +1 \ . S. Hence there is an integer a < 6' such that det(a) 



is non-zero. By equation (5.2 1 we obtain an invertible system of equations 



f t (a)= Y, wt ( a ) ' Na{G, <j>,w). (5-3) 
wew A (G) 



The coefficients N A (G, 4>, w) can be now obtained in polynomial time by Lemma |5.4| This finishes the 
proof of COUNT pin (A) = EVAL pin (yl). The proof for COUNT(A) = EVAL(^) also follows by the 
argument just presented, as the given pinnings remain unaffected. n 

5.2 Dealing with Vertex Weights 

Lemma 5.6 (Theorem 3.2 in [ DGOOD . Let A G IR™ xm be a symmetric matrix with non-negative en- 
tries such that every pair of rows in A is linearly independent. Let D G M™ xm be a diagonal matrix of 
positive vertex weights. Then 

EVAL P '"(A) < EVAL p/n (A D) and EVAL(A) < EVAL(A D). 

After some preparation, we will prove this lemma in Section [5.2.2| 

5.2.1 Some Technical Tools 

The following is a lemma from [DGOO] (Lemma 3.4). 

Lemma 5.7. Let A G IR™ xm be symmetric and non- singular, G = (V,E) a graph, (f> a pinning, and 
FCE.Ifwe know the values 

f r (G) = £ c A (a) J] A r a{u)a{v) (5.4) 

(pCa:V^-[m] uv£F 

2 

for all r G [(\F\ + l) m ], where ca is a function depending on A but not on r. Then we can evaluate 



u)cr(v) 

(j>C(7:V—y[m] uv£F 



in polynomial time. 
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5 EVALUATION AND COUNTING 



Proof. As A is symmetric and non-singular, there is an orthogonal matrix P such that P T AP =: D is 
a diagonal matrix with non-zero diagonal. From A = PDP T we have A r = PD r P T and thus every 
entry of A r satisfies A^ = (PD r P T )ij = Y,™=i p *(v)iM P a{v)v.{ D wY . Hence equation <(5]4j) can be 
rewritten as 



f r (G) = II E P -( 

4>Ca:V—>-[m] uvGF /i=l 



r 



Define the set 

W = I fJ(Ai) a< I < q» for all i € [m], o< = |F| 

U=l i=l 
which can be constructed in polynomial time. We rewrite 



f r (G) = £ * 



for unknown coefficients c w . By interpolation (cf. Lemma 5.4 1, we can recover these coefficients in 
polynomial time and can thus calculate fo(G) = J2 w eW Cw ~ ^ e nave 

m 

MG) = n E p -w^(^(^)° 

4>Ca:V—>[m] uv£F [i=l 

<f>Ccr:V— >-[m] uueF 

which finishes the proof. n 
The following two lemmas are restatements of those in [DGOO] (see Lemma 3.6, 3.7 and Theorem 3.1). 

Lemma 5.8. Let A £ jjmxm ^ e a symmetric matrix in which every pair of distinct rows is linearly 
independent. Let D £ jj mxm ^ e a diagonal matrix of positive vertex weights. Then every pair of rows 
in ADA is linearly independent. Furthermore there is an < e < 1 such that for all i ^ j 



\(ADA)ij\ < eJ (ADA)u(ADA) 



.i.i 



Proof. Define Q = AD^/ 2 \ We have ADA = ADWdW A t = QQ T . That is (ADA)ij = 
(Qi,*, Qj,*} every pair of rows in Q is linearly independent as it is linearly independent in A. Hence, by 
the Cauchy-Schwarz inequality 



(Qi,*-> Qj,*} ^ y (Qi : *i Qi,* 1 ) (Qj,* j Qj,*^) 

which implies that the corresponding 2x2 submatrix of ADA defined by i and j has non-zero determi- 
nant. The existence of e follows. n 

Lemma 5.9. Let A £ ^ mxm }j e a symmetric non-negative matrix in which every pair of distinct rows 
is linearly independent. Let D £ K mxm l, e a diagonal matrix of positive vertex weights. Then there is 
a p £ N such that the matrix 

(ADA)**) 

is non-singular. 



5.2 Dealing with Vertex Weights 
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Proof. Let A' = ADA and consider the determinant 

m 

det(A')= £ sgn^H^) 



i=l 



where S m is the set of permutations of [m]. For some 7r G £ m define t(ir) = \{i \ ir(i) ^ Let e be 
as in Lemma I5T81 Then 



n < £ ,w n ^ n a«».<, = ^ n < 



i=l i=l i=l 

Let id denote the trivial permutation. Then 

/ m \ P 



8=1 



det((A 



> 



vi=l 



By equation (5.5 1, we have 



II 1' 

u=l 



> 



7re5 m \{id} \«=1 / 



7reS m \{id} \i=l / 



(5.5) 



and hence, as < e < 1, the matrix (ADA)^ is non-singular for large enough p 



5.2.2 The proof of Lemma 5.6 



By Lemma 5.9 there is a p G N such that (ADA)^ is non-singular. We will fix such a p for the rest of 
the proof. 

Let G, 4> be an instance of EVAL pin (yl) with G = (V, E) a graph. Construct from G a graph G' as 
follows. Define 

V = {v , ■ ■ ■ | v 6 V,d = d G (v)} 

that is, for each vertex v G V we introduce new vertices. Let ft be the pinning which for every 
v G dom(0) satisfies 4>'(vq) = 4>(v). Let E'" be the set which contains, for each v G V the cycle on 
{^0) • • • , v d G (v)-i}> i- e - {^o^i) • • • j Wd G (u)-i u o} ^ E'" . Let E 1 " be the set of edges, such that each edge 
from E incident with v is connected to exactly one of the v; L . Define E' = E" U E'" and denote by G' p r 
the graph obtained from G' by replacing each edge in E'" with a distinct copy of the graph T PiT to be 
defined next. 

To conveniently define the graph T p>r we will first describe another graph T p . This graph consists 
of two distinguished vertices a and b connected by p many length 2 paths from a to b. Then T PjT is 
the series composition of r many copies of T p . The construction is illustrated in Figure [3] We call a 
and b the "start" and "end" vertex of T p and T p r has start and end vertices induced by the start and end 
vertices of the series composition of T p . 

For a graph H with designated "start" and "end" vertex we denote by ZA,D(i,j', H) the partition 
function Za,d{4>i H) for (ft such that it pins the start vertex of H to i and the end vertex to j. Note that, 
by definition, the vertex weights of the pinned vertices do not occur in this term. 

Claim 1. Let C = (ADA)**), then with X = D^ l l 2 \ we have for all i,j G [m] and r G N, 

Z A!D (i,j;T P;r ) = (XuX^dXCXy)^. (5.6) 
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5 EVALUATION AND COUNTING 



a 




2 A 




Figure 3: The graphs T2 and T2 A . 



Proof. Straightforwardly, 

Z AiD (i,j;T p ) = \y^A ik A kj D kk j = (ADA)$ = C tj 

and therefore 

r r 

ZA,D{i,j;T Pi r) = ^ Y[ Z A , D (a(k),a(k + 1);T P ) ]J D a{k) ^ k) 

<r:[r+l]-*[m] fe=l fc=2 
tr(l)=i, (r(r+l)=j 



<r:[r+l]— t[m] k=l 
tr(l)=i, a(r+l)=j 



k=2 



a:[r+l]->[m] fc=l 
a(l)=i,cr(r+l)=j 



By inspection, the last line equals the right hand side of equation (5.6 1 — as claimed. 



Using the expression just obtained for Z A> i)(i,j; T p r ), we will now rewrite Za,d($, G' p r ). To do so, 
we will, for all v S V, count the indices of the vertices i?o, . . . , Vj G r v \i modulo cLq{v). In particular 
Vd G (v) = vo- Define 7 = U.vedom(4>') D <t>>(v),cj>>{v)- We have for all r G N, 

d G {v)-l 

As the vertices in V are grouped according to the vertices in V we have 

d G (v)-l 

n D <r(v),a(v) = n n d °m,vm- 



Using equation (|5.6|), we further see that for each a : V — > [m], the expression 



d G (v)-l 

n d n n ^w^) 
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turns into 

da(v)-l ((vnv\r\ do(v)—l 



11 D a(Vi)Mvi) ——— — = Y{ 11 {[XCX) ) a{Vi)Mvi+l) 

Thus 

d G (v)-l 

Z A;D ((j)',G' Ptr ) = 7 _1 A^y^) Jl ]J ((XCX) r ) (T („ i ) ff ( ) , i+1 ) 

</>'Co-:V"-s>[m] uu6B" «eV i=0 

= 7 _1 ^ IT 11 ((^C^) r )<j(«)a(t;) 

Given the EVAL(A, D) oracle, we can, in polynomial time, evaluate this expression for every r which 
is polynomial in the size of G. Therefore Lemma [577] implies that we can compute the value 

z = 7 _i n a °{u)v(v) n h{u)o{v) 

4>'Qa:V'^[m]uv£E" uveE'" 

in polynomial time. The proof follows, if we can show that 7 • Z = Za(4>, G). To see this, note that, for 
every configuration a : V -> [m] the corresponding weight n«„eJ5" A v(u)a(v) UuveE'" 4(«)<x(u) in th e 
above expression is zero unless the following holds: For all v 6 V we have a(vo) = . . . = a(vd-i) for 
d = dc{v). for such a configuration, define a configuration a' : V — > [m] such that a(v) = a(vo) for 
every v £ V. It follows from the construction of G' that then 

II A <j(u),a(v) = II A a'{u),<j>{v)- 
uvGE" uvGE 

Since every configuration a' : V — > [m] arises this way, we have 7 • Z = Z a(^i G). This finishes the 
proof of EVAL pm (A) < EVAL pin (A, D). Since the proof is correct also for empty input pinnings 4>, 
this also proves EVAL(A) < EVAL(A, D). 



6 The Pinning Lemma 



In this section we shall give the proof of the Pinning Lemma 4.3 Before we do this, however, we will 
introduce a technical reduction which will also be used later. 



6.1 The Twin Reduction Lemma 

For a symmetric mxm matrix A, we say that two rows Ai^ and Aj t * are twins, if Aj,* = Aj t *. A matrix 
A is twin-free if A4 * 7^ A,,* for all row indices i ^ j. The concept of twins induces an equivalence 
relation on the rows of A. Let 1%, . . . , be the equivalence classes of this relation. The twin resolvent 
of A is the k x k matrix [A] , defined by 

[A]ij = A^ v for some fj, € I, and v £ Ij. 

We say that [A] is obtained from A by fwin reduction. Further, we define the twin resolution mapping 
t : [m] —> [k] of A in such a way that for all \i £ [m] we have /x € ir(p)- That is, r maps every /i € [m] 
to the class L it is contained in. Hence 



[ A ]r(i),r(j) = A i,j for all i,j e [mj. 



(6.1) 
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6 THE PINNING LEMMA 



To use twin reductions in the context of partition functions we need to clarify the effect which twin 
reduction on a matrix A induces on the corresponding diagonal matrix of vertex weights D. We will see 
that the diagonal k x k matrix defined in the following captures this effect 

D\f = A/,„ for all i G [k] . (6.2) 



Lemma 6.1 (Twin Reduction Lemma). Let A e R™ be non-negative and symmetric and D € 
t™ xm a diagonal matrix of positive vertex weights. Let [A] be the twin resolvent of A. Then 



EVAL([A],D^) = EVAL(A,D) and EVAL pin ([A], D^) = EVAL pin (A, D). 

Proof. Let r be the twin-resolution mapping of A. Let G = (V, E) be a graph and <p a pinning. Define 
V' = V \ dom(</>). By the definition of r we have 

Za,d(4>, G) = E II A a{u),c(v) \\ D a(v),a(v) 

<t>Ca:V^[m]uv£E v£V~' 

= E 11 ^roa{u),Toa(v) \\ At(d),ct(u) ■ 

<l>Ccr:V^[m]uv£E veV 

For all configurations a : V — > [m] we have r o a : V — > [k] . Hence, we can partition the configurations 
a into classes according to their images under concatenation with r and obtain 

Za,d(<P,G) = ^ E IT \. A \<r'{u),a'{v) II D c(v),c{v) 

a':V^[k](pCa:V^[m]uveE v£V 

TOO=o' 

e n^'w/w A ( (j ')' 

To4>Ca':V^{k] uv&E 

Here, A(cr / ) is defined by 



a ( (j/ ) = e n D °<?u<?y 

Fix some a' : V — > [k] , we will argue that 



<t>Ca:V^r[m] v€V 
Toa=cr' 



v&V' 



For every configuration a : V — > [m] we have r o a = a' if, and only if, a' 1 ({i}) = a 1 (h) for all 



i G [A;]. Define, for each % G [A;], the set Vi := a' and the mapping fa := <p tdom(</>)nV;- Then 

A ^ = E U d «wmv) 

cf>Ca:V^[m\ v&V 
ViS[fc]: a(Vi)CIi 



6.2 Proof of the Pinning Lemma 
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Define V- := Vi \ dom( ( /) i ), then 



a ^ = n e n^w«) 

i=i vev? ueh 

~ 11 U <j'{v),ct'(v) 



This proves equation (|6.3|). Therefore, 



^A,D(0, G) = ^ ]J [4]a'(u),cr» II D £(v),a>(v) = Z [A],DW ( r ° & G ). 

TO(j>C(T>:V^[k] uv&E vGV 

This witnesses the claimed reducibilities. 
6.2 Proof of the Pinning Lemma 



We shall prove the Pinning Lemma 4.3 now. The proof of this Lemma relies on a result of Lovasz 
BLov0611 . To state this result we need some further preparation. Let A, D be m x m matrices and A', D' 
be n x n matrices such that D and D' are diagonal. The pairs (^4, D) and (^4', D 1 ) are isomorphic, if 
there is abijection a : [m] ->• [n] such that Aij = ^,(i)a(j) for a11 G t m ] and ^ = ^(i),^) for 
all i € [m]. An automorphism is an isomorphism of (y4, Z?) with itself. 

We will moreover need to consider pinnings a bit differently than before. Rather than defining a 
pinning 4> for some given graph G, it will be convenient in the following to fix pinnings and consider 
graphs which are compatible with these. To define this more formally, we fix tfi : [k] — > [m] to denote 
our pinning. A k-labeled graph G = (V, E) is then a graph whose vertex set satisfies V D [k]. In this 
way (j> is compatible with every /c-labeled graph. 

Lemma 6.2 (Lemma 2.4 in HLov061 ). Let A G K mxm be a non-negative symmetric and twin-free ma- 
trix and D G K mxm a diagonal matrix of positive vertex weights. Let </>, ip be pinnings. If Za,d{4>-, G) = 
Za^d^i G) for all k-labeled graphs G, then there is an automorphism a of (A, D) such that eft = a{ip). 

Utilizing this result, we can now prove the pinning result first for twin-free matrices. 

Lemma 6.3. Let A G M™ xm be non-negative, symmetric, and twin-free and D G M™ xm a diagonal 
matrix of positive vertex weights. Then 

ENAL pin (A, D) = EVAL(A D). 

Proof. As EVAL(A, D) < EVALP' n (yl, D) holds trivially we only need to prove EVAL pin (A, D) < 
EVAL(A, D). 

Let G = (V, E) and a pinning (f> be an instance of EVAL pm (^4, D). By appropriate permutation 



of the rows/columns of A and D (cf. Lemma 5.3) we may assume that [k] = img (f> C [m] for some 



k < m. Let G = (V, E) be the graph obtained from G by collapsing the the sets for all i e [k]. 

Formally, define a map 

i ,d£ for some i € [k] 

v , otherwise 



70) 
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6 THE PINNING LEMMA 



Then G is a /c-labeled multigraph (with possibly some self-loops) defined by 

V = [k] U (V \ dom(0)) 
E = {7(^)7(7;) I uv G E}. 

Recall that in partition functions of the form Za,d{(I>, G) vertices pinned by <fi do not contribute any 
vertex weights. Hence, Za,d(4>, G) = (idrw , G) where idny denotes the identity map on [k]. Call 
two mappings x, ip ■ [k] ^ [m] equivalent if there is an automorphism a of (A, D) such that x = ct o tp. 
Partition the mappings : [k] — > [m] into equivalence classes 1% , . . . , I c according to this definition and 
for all i G [c] fix some tpi G Ij. Assume furthermore, that tpi = idny. Clearly for any two x, ip from the 
same equivalence class, we have Za,d(x> F) = ■^ j 4,d(V ; ) F) f° r every graph F. Therefore, for every 
graph G', 

Za,d(G') = ZaM^u G') ■ I J2 II I ( 6 - 4 ) 

i=l V'/'G-fi ti£dom(i/)) / 



Define, for each i G [c] the value q = (^E^e/, Ilvedom^) D ^{v)ti>{v)) ■ We claim the following 

Claim 1. Let / C [c] be a set of cardinality at least 2 such that 1 € I. Assume that we can, for every 
/c-labeled graph G' , compute the value 

Y,Ci-Z A ,D{ipi,G'). (6.5) 

Then there is a proper subset I' C I which contains 1 such that we can compute, for every /c-labeled 
graph G" , the value 

Y^^-ZaM^G"). (6.6) 

This claim will allow us to finish the proof. To see this, note first that by equation ( |6.4| ) we can compute 
the value ( |6.5j ) for I = [c] and G' = G. Thus after at most c iterations of Claim [I] we arrive at 
ci • Za,d (tpi , G) . Further, c\ is effectively computable in time depending only on D and therefore we can 
compute Z A ,D(id [k] , G) = Z a ,d($, G). This proves the reducibility EVAL pin (A, D) < EVAL(A, D). 



Proof Of Claim[l| Assume that we can compute the value given in (6^ I. Lemma [672] implies that for 
every pair i ^ j £ I there is a /c-labeled graph F such that 

Z A Mi,T)^ZA,D{^,T). (6.7) 

Fix such a pair i 7^ j G I and a graph T satisfying this equation. Note that this graph can be computed 
effectively in time depending only on A, D and fy, ipj. Let G s denote the graph obtained from G by 
iterating s times the /c-labeled product of G with itself. We can thus compute 



^ ■ z AD (i>,, GT-) = ]T c,z A Mi>>, a') ■ z A:D y,i,Ty. (6.8) 



E 

iei iei 



Partition I into classes Jq, . . . , J z such that for every v G [0, z] we have i',f G J v if, and only if, 
Za,d (ipi 1 > r) = ZA,D(ipj')T). Since one of these sets J v contains 1 and all of these are proper subsets 
of /, it remains to show that we can compute, for each v G [0, z], the value 

ciZ A ,D(ipi>,G'). 
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To prove this, define x u := Za,d{^v, T) for each v £ [z] and an i' G J y . Equation ( |6.8| ) implies that 
we can compute 

One of the values might be zero. Assume therefore w.l.o.g. that xo = 0, then evaluating the above 
for s = 1, . . . , z yields a system of linear equations, which by Lemma 5.4 can be solved in polynomial 
time such that we can recover the values J2i'eJ v c i' z A,D{^v ■, G') for each v > 1. Using equation ( |6.5| ) 
we can thus also compute the value 

I Z~2, Ci ' Z A,D(lpi,G') I ~E E c i' Z A,D{^i',G') = ^ c i ,Z A,D($i'iG'). 

\iei J v=\ \i'ej v J i'eJo n 



The proof of the Pinning Lemma now follows easily from Lemmas |6.1| and 6.3 



Proof (of the Pinning Lemma 4.3 1. Fix the twin resolvent A' = [A] of A and let D' = D^ A \ It suffices 
to show 

EVAL pin (A', D') = EVAL(,4 / , D'). (6.9) 
To see this note that by the Twin Reduction Lemma [6TT] we then have the chain of reductions 

EVAL pin (A, D) = EVAL pin (^4 / , D') = EVAL(A', D') = EVAL(A, D). 

The proof of equation \6.9) follows from Lemma [63] n 



7 The General Conditioning Lemma 
7.1 Dealing with {0, 1} Matrices 

As a technical prerequisite, we need a part of the result of [DG00] (Theorem 1.1 in there). We call a 
block non-trivial if it contains a non-zero entry. 

Lemma 7.1 (#H-Coloring Lemma). Let Abe a symmetric connected and bipartite {0, l}-matrix with 
underlying non-trivial block B. If B contains a zero entry then the problem EVAL f ""(A) is j^V-hard. 

Proof. We will start with the following claim which captures the main reduction. Call a matrix A 
powerful if it is a symmetric connected and bipartite {0, 1} -matrix whose underlying block B contains 
a zero entry. 

Claim 1. Let A be a powerful matrix with underlying m x n block B. If either n > 2 or m > 2 then 
there is a powerful matrix A' with underlying m! x n! block B' such that 2 < m' < m and 2 < n' < n, 
at least one of these inequalities is strict and 

EVAL pin (yl / ) < EVAL pin ( J 4). 

Before we prove the claim, let us see how it helps in proving the lemma. Let A be as in the statement 
of the lemma. Iterating Claim[T]for a finite number of steps we arrive at EVAL pin (A / ) < EVAL pin (^4) 
such that the block B' underlying A' is a 2 x 2 matrix of the form (up to permutation of rows/columns) 
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For every empty pinning <j) and bipartite graph G with c connected components Za'{4>, G) equals 2 C 
times the number of independent sets of G. This straightforwardly gives rise to a reduction from 
the problem of counting independent sets in bipartite graphs which is well-known to be #P-hard (see 
[PB81). 



Proof of Claim[I] As B is a block with a zero entry, there are indices i, j, k, I such that B^ = Bu = 
Bjk = 1 and Bji = 0. Fix these indices and let I = {u \ Bi v = 1} and J = {fj, \ B^ = 1} be the sets 
of indices of 1-entries in row i and column k. Let A* be the connected bipartite matrix with underlying 
block Bjj. We will show that 

EVALP in (,4*) <EVALP in (,4). (7.1) 

To see this, let G, be an instance of EVALP'"^*). We shall consider only the case here that G = 
(U U W, E) is connected and bipartite with bipartition U, W and that <j) P ms a vertex a G U to a row of 
Bjj. All other cases follow similarly. 

Define G' as the graph obtained from G by adding two new vertices v! and w' and connect every 
vertex of U by an edge to w' and every vertex of W by a edge to v! . Let 4>' be the pinning obtained from 
<fi by adding u' i->- i and w' h-> fc. We have Za* (4>, G) = Za{4>' i G') which yields a reduction witnessing 
equation ( |7.1| ). 

Before we proceed, we need another 

Claim 2. Let A + be a powerful matrix with underlying m + x n + block £> + . There is a twin-free 
powerful matrix A" with underlying m" x n" block S" such that 2 < m" < m + and 2 < n" < n + and 

EVAL^^") < EVAL pin (^ + ). 



Proof. This is a straightforward combination of the Twin Reduction Lemma 6.1 and Lemma 5.6 



Combining equation ( |7T) and Claim |] we arrive at EVAL pin (A // ) < EVAL pin (A) for a powerful 
twin-free matrix A". The block B" underlying A" has some dimension m" x n" such that 2 < m" < m 
and 2 < n" < n. Further, up to permutation of rows/columns, this block has the following form 



/ 1 



B" 



1 \ 



V i * ... / 

To prove Claim[T]it suffices to devise a reduction witnessing 

EVAL pin (^) < EVALP 1 ^^") 



(7.2) 



for a powerful matrix A' such that the block B' underlying A' has either fewer rows or columns than 
B". 

To devise such a reduction, assume that B" has at least three rows (the case that B" has at least 
three columns is symmetric). As B" is twin-free, there must be entries B'^ a 7^ B^ a . Assume w.l.o.g. 
that -Bga = an d = !• ^ s B" i s twin-free we further have that B'^ b = for some b. Let 
K = {k I i?3 K = 1} be the indices of non-zero entries of B'^. Define B' = B" K and let yl' be the 
connected bipartite matrix with underlying block B'. It remains to devise the desired reduction. 

Let G,4> be an instance of EVAL pin (A / ). As before, we shall consider only the case that G = 
(U U W, E) is connected and bipartite with bipartition U, W and that <j) P ms a vertex a G £7 to a row of 
5'. All other cases follow similarly. 



7.2 From General Matrices to Positive Matrices 
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Define G' as the graph obtained from G by adding one new vertex u' and connect every vertex of 
W by a edge to u'. Let <fi' be the pinning obtained from <fi by adding u' i— > 3. This yields a reduction 
witnessing equation ( |7.2[ ) and hence finishes the proof of Claim [T] n 



7.2 From General Matrices to Positive Matrices 



In the first step of the proof of the General Conditioning Lemma 4.4 we will see how to restrict attention 
to matrices with positive entries. 

Lemma 7.2 (The Lemma of the Positive Witness). Let A G E> mxm be a symmetric non-negative ma- 
trix containing a block of row rank at least 2. Then there is an ^-matrix A' satisfying condition (A), that 
is A' is positive symmetric of row rank at least 2, such that 

EVAL^'V) < EVAL p/n (,4). 

The proof of this lemma relies on the elimination of zero entries. We do this in two steps, first (in the 
next Lemma) we pin to components. Afterwards, we will see how to eliminate zero entries within a 
component. 

Lemma 7.3 (Component Pinning). Let A 6 g mxm t, e a symmetric matrix. Then for each component 
C of Awe have 

EVAL p/n (C) < EVAL p/n (,4). 

Proof. Let G,<p be the input to EVAL pin (C) for some component C of some order m! x ml. Let 
Gi, . . . , Gi be the components of G and 4>i, - ■ ■ ,<Pe tne corresponding restrictions of <fi. Since 



Z c (<t>,G) = f[Z c (cf>i,G i 



we may assume w.l.o.g that G is connected which we shall do in the following. Let v be a vertex in G 
which is not pinned by (f>. Define, for each i 6 [m], <fii as the pinning obtained from <p by adding v i— > i. 
We have 

m 

Z c (<!>,G) = J2 z c(<l>i,G). 

i=l 

Thus it suffices to compute Zc{4>i, G) which can be obtained straightforwardly as these values equal 

ZA{(f>i,G). □ 

Lemma 7.4 (Zero-Free Block Lemma). Let A <E E> mxm be a symmetric matrix. Then either EVAL pin (A) 
is j^P-hard or no block of A contains zero entries. 

Proof. Let A" be the matrix obtained from A by replacing each non-zero entry by 1. Let C be a 
component of A" whose underlying block B contains a zero entry. For every given graph G and pinning 
4> each configuration a which contributes a non-zero weight to Z a(<P, G) contributes a weight 1 to 



Z A/l ((j),G). Thus by EVAL pin (A) = COUNT pm (^l) (cf. Lemma \5.5) and component pinning (cf. 
Lemma [73] > it is easy to see that 

EVAL pin (C) < EVAL^V') < EVAL pin (A). 



If C is bipartite, the result follows from the #H-Coloring Lemma [TT] 
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If C is not bipartite, we need a bit of additional work. Note that in this case C = B. Let A 1 be the 
connected bipartite 2r x 2r matrix with underlying block B. If we can show that 

EYALP m (A') < EVAL pin (£) 

then the result follows from the #H-Coloring Lemma [7TT] 

Let G, <f) be an instance of EVAL pin (A / ). For simplicity, assume that G = (U U W, E) is connected 
and bipartite. If 4> is empty, then we have Z A /(G) = 2 • Zb(G). Otherwise, Z A i(<j), G) = unless 4> 
maps all elements of dom(0) n U to [r] and all entries of dom(0) n W to [r + 1, 2r] or vice versa. Since 
both cases are symmetric, we consider only the first one. Let <p' : dom(</>) — > [2r] be the pinning which 
agrees with <p on dom(0) n U and for all w G dom(0) n satisfies <^'(u;) = </»(iu) — r. By inspection 
we have 

Z A ,( ( f>,G) = Z B ^',G'). n 



Proof ( of Le mma [Z2| ). Let A G § mxm contain a block of ran k at least 2. By 2-stretching (cf. 
Lemma 5.1 1 we have EVAL pin (^4 2 ) < EVAL pm ( J 4) and Lemma 5.2 gu aran tees that A 2 contains a com- 
ponent BB T which has rank at least 2. By component pinning (Lemma 7.3 1 we have EVAL pm (BB T ) < 
EVAL pin (yl 2 ). If BB T contains no zero entries, we let A' = BB T . 

Otherwise let A ' G N 2x 2 be a matrix satisfying the conditions of the lemma. As BB T contains a zero 
entry, Lemma flA\ implies that EVAL pin ( J B5 T ) is #P-hard and thus EVAL pin (A') < E\AL^ n (BB T ) 



finishing the proof. n 
7.3 From Positive Matrices to X-matrices 

Let A be a matrix which satisfies condition (A). We will now see, how to obtain a matrix A' which 
additionally satisfies (B). That is, we will prove the following lemma. 

Lemma 7.5 (X-Lemma). Let A G § mxm be a matrix satisfying condition (A). Then there is an S- 
matrix A' satisfying conditions (A) and (B) such that 

ENAL pin {A') < EVAL p/n (,4). 

To prepare the proof, we present some smaller lemmas which will also be useful in later sections. 

Lemma 7.6 (Prime Elimination Lemma). Let A G Z mxm (A e Z[X] mxm ) and pbe a prime number 
(an irreducible polynomial). Let A' be the matrix obtained from A by replacing all entries divisible by 
p with 0. Then 

EVAL pin (A') < EVAL pin (A). 



Proof By Lemma |5.5| it suffices to give a reduction witnessing EVAL pm (yl / ) < COUNT pm (A), which 
we will construct in the following. Let G = (V,E) be a given graph cj) a pinning. As Wa'{G) = 
Wa(G) \{w I w divisible by p}, we have W A '(G) C W A (G). Moreover N A >(G, (f), w) = N A (G, (f), w) 
for all w G W A > (G) which implies 

Z A ,(<j>,G)= J2 w-N A {G,<j>,w). 
wew A ,(G) 

The values N A (G, 4>, w) can be obtained directly using the COUNT pm ( J 4) oracle. n 



7.3 From Positive Matrices to X-matrices 
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Let p be a prime number (an irreducible polynomial, respectively) and a G Z (a G Z[X], resp.)- Define 

^max{fc>0|p fe divides a} jf q ^ Q 

, otherwise. 

For a matrix A the matrix A\ p is then defined by replacing each entry Aij with Aij\ p . 

Lemma 7.7 (Prime Filter Lemma). Let A G Z mxm (A G Z[X] mxm ) and p be a prime number (an 
irreducible polynomial). Then 

EVXL pin (A\ p ) < EVXL pin (A). 



Proof. By Lemma 5.5 it suffices to give a reduction witnessing EVAL pin (A| p ) < COUNT pin (,4), which 
we will construct in the following. For a given graph G = (V, E) and a pinning <p we have 

Z A \ P (<P,G)= w\p-Na{GA,w)- 

wdW A {G) 

The values Na(G, 4>, w) can be obtained directly using a COUNT pin (^4) oracle. n 

Lemma 7.8 (Renaming Lemma). Let p G Z[X] \ {-1,0,1} and A G Z[X] mxm a p-matrix. Let 
q G Z[X] and define A' G Z[X] mxm by 



j^i . _ ( q l , f/iere ij an I > s.t. A^ = p l 
% i \ , otherwise 



That is, A' is the matrix obtained from A by substituting powers of p with the corresponding powers of 
q. Then 

EVAL pin (A') < EVAL pin (A). 

Proof. Consider A' as a function A'(Y) in some indeterminate Y. We have A = A'(p). Let £ max be 
the maximum power of Y occurring in A'(Y). For every graph G and pinning <p, the value f(Y) := 
Z A HY){<fii G) is a polynomial in Y of maximum degree \E\ ■ £ max . By i-thickening (cf. Lemma 5.1 ), 
using an EVAL pm (yl) oracle, we can compute the values f(p t ) for t = 1, . . . , \E\ ■ £ max - Thus, by 
interpolation (cf. Lemma 5.4) we can compute the value f(q) for a q as given in the statement of the 
lemma. By f(q) = Z^i^((j), G) this proves the claimed reducibility. n 

Lemma 7.9 (Prime Rank Lemma). Let A G Z mxm (A G Z[X] mxm ) contain a block of row rank at 
least 2. There is a prime number (an irreducible polynomial) p such that A\ p contains a block of row 
rank at least 2. 

Proof. Let A^* and Ay^ be linearly independent rows from a block of A. Assume, for contradiction, 
that for all primes (irreducible polynomials, resp.) p every block in A\ p has rank at most 1. We have, for 
all j G [m] , 

Aij = J_ J_ \p and A^j = A^j \ p 
p p 
where the products are over all primes (irreducible polynomials, resp.) dividing an entry of A. By 
assumption, there are a p , (3 p G Z (in Z[X], resp.) such that a p ■ Ai *\ p = j3 p • -Aj'*|p for all primes 
(irreducible polynomials). Therefore, 

Aj n a p = n a P Ai i \p = n /w? = ^ ■ n ^ f ° r &n j g h ■ 

p p p p 

And hence, Ay * and Aj * are linearly dependent — a contradiction. n 
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Dealing with algebraic numbers. The following lemma tells us that the structure of the numbers 
involved in computations of partition functions on matrices with algebraic entries is already captured by 
matrices with natural numbers as entries. 

Lemma 7.10 (The Arithmetical Structure Lemma). Let A G M™ xm be symmetric and non-negative. 
There is a matrix A' whose entries are natural numbers such that 

EVALP in (A') = EVALP in (A). 

If further A contains a block of rank at least 2 then this is also true for A'. 

We have to introduce some terminology. Let B = . . . , b n } C M A be a set of positive numbers. The 
set B is called multiplicatively independent, if for all Ai, . . . , A n G Z the following holds: if 6 Al • ■ ■ 6 An 
is a root of unity then Ai = . . . = A n = 0. In all other cases we say that B is multiplicatively dependent. 
We say that a set S of positive numbers is effectively representable in terms of B, if for given x G S we 
can compute Ai, . . . , A n G Z such that x ■ 6 Al ■ ■ ■ & An = 1. A set B is an effective representation system 
for a set S, if S is effectively representable in terms of B and B is multiplicatively independent. 
We need a result from [RicOl ] which we rephrase a bit for our purposes. 

Lemma 7.11 (Theorem 2 in I RicOll ). Let ai, . . . , a n G Q(9) be positive real numbers given in stan- 
dard representation, each of description length at most s. There is a matrix A G Z nx ™ such that, for 
vectors A G Z n we have 

n 

Y[ a Al =lif, and only if, A ■ A = 0. (7.3) 

i=l 

The description length of A is bounded by a computable function in n and s. 

This result straightforwardly extends to an algorithm solving the multiplicative independence problem 
for algebraic numbers. 

Corollary 7.12. Let a%, . . . , a n G Q(9) be positive reals given in standard representation. There is an 
algorithm which decides if there is a non-zero vector A = (Ai, . . . , A n ) G Z n such that 

n 

U a i i = 1 ( 7 - 4 ) 

i=i 

Furthermore, if it exists, the algorithm computes such a vector A. 

Lemma 7.13. Let S C M4 be a set of positive numbers. There is an effective representation system 
B C IR^ of positive numbers for S which can be computed effectively from S. 

Proof. We shall start with the following 

Claim 1. If S is multiplicatively dependent then there is a set B' C M A of non-negative numbers such 
that \B'\ < \S\ and S is effectively representable by B' . 



such that b^ 1 ■ ■ ■ b^ n = 1. We can easily make sure that at least one of the Aj is larger than zero. Assume 
therefore w.l.o.g. that Ai > 0. Fix a set B' = {b' 2 , ■ ■ ■ , b' n } where each b' { is the positive real Ai-th root 
of bu that is (6;) Al = bi. Then 



Proof. Let S = {61, ... , b n } then Corollary 7. 12 implies that we can compute a non-zero vector A G 



1/ 

Ai 



&1 1 • Yl( b 'i) Xl = 1 and hence h ■ J|(6^) Al = 1. 

\i=2 J i=2 

All operations are computable and effective representation of 5 by B' follows. 



7.3 From Positive Matrices to X-matrices 
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Apply Claim [T] recursively on S. Since the empty set is multiplicatively independent, after at most 
finitely many steps, we find an effective representation system B for S. n 



Proof (of Lemma 7.10 ). Let S be the set of non-zero entries of A. By Lemma 7.13 we can compute 
an effective representation system B for S. However, with respect to our model of computation we 
need to be a bit careful, here: assume that S C Q(6>) for some primitive element 9. The application 



of Lemma 7.13 does not allow us to stipulate that B C Q(6). But in another step of pre-computation, 
we can compute another primitive element 0' for the elements of B such that B C Q(6') (c.f BCoh 93P. 
Then we may consider all computations as taking place in Q(9'). 

Assume that B = {b±, . . . , &„}, then every non-zero entry of A has a unique computable represen- 
tation 



V 



Let pi, . . . ,p/3 be j3 = \B\ distinct prime numbers and define A' as the matrix obtained from A by 
replacing in each non-zero entry Aij the powers of b G B by the corresponding powers of primes, that 
is, 



A'ij \[i> 



Recall the definition of Wa{G) in equation (5J_ ). For each w G Wa(G) we can, in polynomial time 
compute a representation w = Yiij A ij J as powers of elements in S. The effective representation of S 
in terms of B extends to Wa(G) being effectively representable by B. Moreover, as S depends only on 
A, the representation of each w G Wa(G) is even polynomial time computable. We have 

Z A {<t>,G)= w-N A (G,^w) 

w&A) A (G) 

In particular, for each w G Wa(G), we can compute unique A TO) i, . . . , X w ,n G Z such that W'b-^ ' • • • bn 
1. Define functions / and g such that for every w G Wa(G) we have 

n n 

f(w) = l[pi-" and g( w ) = l\b X -. 

v=\ v=l 



Thus we obtain 

/H 

This yields a reduction for EVAL pin (yl / ) < EVAL pin (,4). The other direction follows by 

/H 



Za>(<P,G)= V w^-N A (G,<p,w). 



z A (<P,G)= w ' -jrl-NAGA^'). 



This also proves the reducibilities EVAL(A') = EVAL(yl) since the input pinnings remain unaffected. 
To finish the proof it remains to consider the case that A contains a block of rank at least 2. We have 
to show that A' has this property as well. Let us now argue that A' contains a block of rank at least 
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2. Let Ai t * and A41 * be linearly independent rows from a block of A. Assume, for contradiction, that 
A'-, „ = a • A'- for some a. Let J be the set of indices j such that A'- - ^ 0. For each 7 G J we have 



Hence, for (3 



a 



A V j ■ (A\ 



13 > 



n 



b X n' 3n Xl3n we obtain A t 



ft ■ Ai^ 



a contradiction. 



7.3.1 Proof of the X-Lemma|73J 

Let A £ § mxm be a matrix satisfying condition (A). Recall that S is one of R A , Q, Z, Z[X] or ( 

If S = Ma then the entries of A are all positive real values. Thus Lemma 7.10| implies that there is a 

positive matrix A' E Z mxm of rank at least 2 such that EVALP' n (yl / ) = EVALP' n (,4). 

If A is a matrix of entries in Q (Q[X], respectively) then let A be the lowest common denominator 

of (coefficients of) entries in A. For a given graph G = (V, E) and pinning <p we have 

Z XA (^G) = \\ E \Z A (^G). 

The mattix A • A is a matrix with entries in Z (Z[X], respectively). 

It remains to prove the Lemma for the case that S is either Z or Z [X] . By the Prime Rank Lemma 7.9 
there is a prime (irreducible polynomial) p such that A\ p contains a block of rank at least 2. Fix such a 
p and define A' = A\ p . By the Prime Filter Lemma 



7.7 



we have 



EVAL^V) < EVAL pin (j4). 



Furthermore, by the Renaming Lemma 7.8 we may assume that A' is an X-matrix. This completes the 
proof. 



7.4 From X-matrices to the General Conditioning Lemma 

Lemma 7.14. Let A £ Z[X] mxm be symmetric and positive such that not all 1-entries of A are con- 
tained in 1-cells. Then EVAL pm (A) is j^V-hard. 

Proof. If not all 1-entries are contained in 1-cells, then there are i, j, k, I such that A^ = An = A^ = 1 
and Aji ^ 1. Let A' be obtained from A by replacing eac h en try not equal to 1 by 0. We have 



7.6 



EVALP |n (^') < EVALP |n (^) by the Prime Elimination Lemma[ 

By construction A^ = An = A jh = 1 and Aji = 0. Thus A' contains a block with zero entries and 
EVAL pin (A') is #P-hard by Lemma" 



7.4 



Proof (Of the General Conditioning Lemma \4~4\ . By Lemma 7.2 ther e is a matrix C' which satisfies 
(A) such that EVALP' n (C") < EVAL pin (^4). Then by the X-Lemma 
satisfies (A) and (B) such that EVALP' n (C) < EVALP' n (C"). Lemma" 



7.5 



7.14 



there is a matrix C which 
implies that EVALP' n (C) is 
#P-hard if not all 1-entries are contained in 1-cells. In this case we let A 1 be some 2x2 matrix satisfying 
conditions (A)-(C). 

Assume therefore that all 1-entries of C are contained in 1-cells. If C contains exactly one 1-cell 
then the proof follows by the symmetry of C. We only have to make sure that we can permute the 



entries of C such that condition (C) is satisfied. This is guaranteed by the Permutability Principle [53 
Assume therefore that C contains more than one 1-cell. Define an X-matrix C 



C 2 \ 



x- 
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Claim 1. For every 1-cell Ckl of C the principal submatrix C* KK is a 1-cell of C*. 

Proof. Note that C 4 * = 1 only if there is an I such that Cu = Cj£ = 1. This proves the claim. h 

Now C* has all 1-entries in principal 1-cells and is of rank at least 2 by the fact that there are at least two 
1-cells in C. A reduction witnessing EVAL pm (C*) < EVAL pin (C) is given by applying 2-stretching 



(cf. Lemma 5. 1 1 and the Prime Filter Lemma [7/7] in this order. n 



8 The Two 1-Cell Lemma 

The (TIC) - Conditions for Matrices with two 1-cells. We define two additional conditions. 
(TIC -A) A has at least two 1-cells. 

(TIC - B) All diagonal entries of A are 1. 

A l-row (l-column) in a matrix A is a row (column) which contains a at least one 1 entry. We call all 
other rows (columns) non-l-rows (non-l-columns). 

Lemma 8.1 (1-Row-Column Lemma). Let A G ^x] mxm be a positive and symmetric X-matrix. Let 
A' be obtained from A by removing all non-l-rows and non-l-columns. Then 

EVAL p/n (y4 / ) < EVAL p/n (^). 

Proof. For every i G [m] let Cj denote the number of 1-entries in row A^. And let A" be the matrix 
defined by A'[- = CiCjA'- for all i, j G [to]. 

We start with a reduction witnessing EVAL pin ( J 4") < EVAL pin (,4). Let G = (V,E), cf) be an 
instance of EVAL pin (yl"). Let A be the maximum degree of X in A, let k = A • \E\ + 1 and define a 
graph G' = (V',E')by 

V' = {v, v 1 \veV} 

E' = £U{ei,...,e* | v G V, e\ = vv 1 , Vt G [k]}. 



We have 



0Co-:V'->[m] uvEE' 

= ^ IT A <?H,a(v) JJ( J 4 CT (t ) ), CT (t) 1 )) fc 

4>C<T:V'^[m] uv£E vGV 

That is, for every a the degree (as a polynomial in X) of the weight expression 

uv£E v£V 

is smaller than k if, and only if, A cr ( v \ cr r v i\ = 1 holds for every v G V. That is, if a configuration a maps 
a vertex of V to a non- l-row, then the degree of the weight expression of a will always be at least k. On 
the other hand, for each configuration a mapping G to A there are exactly ^ ■ ...ct many 
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configurations a' : V' — > [m] of G' extending a in such a way that their weight is of degree smaller than 
k. 

By the polynomial time equivalence of COUNT pin (^) and EVAL pin (A) it suffices to reduce EVAL pin (yl") < 
COUNT pin (A). We can do this by computing the values w ■ Na(G', 4>, w) for all weights w of degree 
smaller than k. Then, with A' = A"\ x , the remaining step EVAL pin (A') < EVAL pin (^4") is given by 
the Prime Filter Lemma [777] n 



Lemma 8.2. Let A G TL\X\ nxm be a matrix which satisfies (A) and contains at least two 1-cells. Then 
there is a matrix A' satisfying conditions (A)-(C) and both (TIC) conditions such that 

EVALP in {A') < EVALP in (A). 

Proof. Let i,j,i',f G [m] be witnesses for the existence of two 1-cells in A in the sense that Aij = 
Ai'ji = 1 but Aj/A ^ 1 and let p be an irreducible polynomial which divides A^j. Let B be the 
matrix obtained fro m A\ v by replacing all pow ers o f p with the corresponding powers of X. Then by 
Prime Filter Lemma ffj\ and Renaming Lemma |7j} we have EVAL pin ( J B) < EVAL pin (A). Note that B 



satisfies condition (B) and it satisfies (A) as A does. We may assume that all 1 -e ntries of B are contained 
in 1-cells, because otherwise EVAL pin (B) would be #P-hard by Lemma ' 



7.14 



2 



\X- 



As it could possibly be the case that not all 1-cells of B are on the diagonal, we form B' = B 
Claim 1. For every 1-cell Bjcl of B the principal submatrix B' KK is a 1-cell of B' . 
Proof. Note that B'^ = 1 only if there is an I such that Bn = Bj£ = 1. This proves the claim. _| 

By this claim, since B contains at least two 1-cells, the matrix B' does so as well and therefore it satisfies 
(A). Condition (B) is satisfied by definition. We have EVAL pin (5') < EVAL pin (B) by application of 2- 



stretchi ng (c f . Lemma pT| ) and the Prime Filter Lemma [777] in this order. Applying the 1 -Row-Column 
Lemma 



on B' then yields EVAL pm (y4') < EVAL pin ( J B / ) for a matrix A 1 which (up to permuting 
rows and columns) has the desired properties. n 

The cells of a matrix A satisfying conditions (A)-(C) are the submatrices Aj i j j for Ij, Ij as defined in 
condition (C). We call such a matrix A a cell matrix, if in each of its cells Aj.j. all entries are equal. 

Lemma 8.3. Let A G 7 J [X] mxm satisfy conditions (A) - (C) and both (TIC) conditions. Then there is 
a cell matrix C G 7 J [X] mxm which also satisfies (A) - (C) and both (TIC) conditions such that 

EVAL p/n (C) < EVAL p/n (^). 

Proof. We define a sequence of matrices with Aq = A and for all v G N we let A v +\ = A 2 v \x- As A 
satisfies conditions (A) - (C) and both (TIC) conditions, this is true for all matrices in the sequ ence . 
Further we have EVAL pin (^4y + i) < EVAL pin (^4y) for all v, by applying 2-stretching (cf. Lemma 
and the Prime Filter Lemma 1777] in this order. 



5.1l 



By definition, deg((^4. I/+ i) i j) = min {&eg{{A v ) ik ) + deg((A u ) jk ) \ k G [to]} foralH, j G [to]. As 
condition (TIC - B) implies deg(A,j) = for all j G [to], we obtain 

deg((A, + i)#) < deg((A u )ij) for all v. 

That is, the degrees of the entries of A v are non-increasing with v and thus there is a \i such that 
A^+x = Aft. To finish the proof we will show that C = A^ is a cell matrix. Let Cjj be a cell of C 



8.1 #P '-hardness 
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which is not a 1-cell. Let i £ I and j G J be such that deg(Cjj) is minimal. By definition we have 
C = C 2 \x and therefore, for all j' G J, 



deg(CV) = deg((C 2 | x ) 



min{deg(Q fc ) +deg(C jk ) \ k G [m]} < deg(Cjj). 



The inequality follows from the fact that Cjj = 1. Then by the minimality of deg(Cij) we have 
deg(CV,/) = deg(Cjj). Analogous reasoning on i' G I yields deg(CV,) = deg(CVj/) for all i' G I, j' G 
J. Thus C is a cell matrix. n 



8.1 #P-hardness 

Lemma 8.4. Let 5 G N araJ A' be a matrix of the form 

( 1 ••• 1 2 5 



A! 



1 

2* 



1 2 d 



\ 2 s 



■ 2 s \ 

■ 2 s 

■ 1 

• 1 / 



Then EVAL p/n (,4) is #P-hard. 



Proof. Let [A'] be the twin resolvent of A'. By the Twin Reduction Lemma 6. 1 we obtain EVAL( [A'] , D) 
EVAL(A') where, for a,b,5 > 1, the matrices [A'] and D satisfy 



[A>] 



1 2 d 
2 s 1 



and!? 



a 
b 



We have EVAL([^4']) < EVAL([A'],D) by Lemma [5^6| Therefore it remains to show that EVALQA']) 
is #P-hard. To see this, let G = (V, E) be a graph and for two sets of vertices U,W cy let e(U, W) 
denote the number of edges in G between U and W. We have 



£ n ^w^w 

er:V->-[m| Uue-E 



Z A (G) 

Let Q be the number of a : V — >• [to] with weight 2 d ' then we have. 

Z A (G) 



^ 2 5-e(<x- 1 (l)^- 1 (2)) 



\E\ 
j=0 



By EVAL(^4) = COUNT(A) we can determine the coefficients c,;. Let u be maximum such that c v / 0. 
Then c v is the number of maximum cardinality cuts in G. Therefore, this yields a reduction from the 
problem #MAXCUT of counting maximum cardinality cuts — a problem well known to be #P-hard 
(this follows, for example, from the work of Simon MSim77IO . n 

Lemma 8.5. Let A G 7 J [X] mxm he a cell matrix satisfying conditions (A)-(C) and both (TIC) condi- 
tions. Then FNAL pin (A) is #P-hard. 
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Proof. Let 5 = minjdeg^y) \ Aij / 1, i, j G [m]} and A = max{deg(^4jj) | Aij / 1, i, j G [m]}. 
Let Aij be a cell of A with entries X s and define A' = ^4(/uJ)(7uJ) which, by symmetry, the definition 
of the cells and (C) has the form 



/ 1 



A' = A 



(/UJ)(/UJ) 



1 

X s 



\ X 6 



1 X s 



1 X s 
X s 1 



X 6 



X s \ 

X s 

1 

1 J 



We will show first that EVALP' n (^') < EVAL pin (A). For a graph G = (V,E) and a pinning <j> of 
EVAL pin (v4'), define a graph G' = (V, E') as follows. Let k = \E\ ■ A + 1, we obtain G' from G by 
adding two apices which are connected to each vertex of V by edges with multiplicity k, that is, 

V' = VU{x,y} 

E' = E U {{xv) k , (yv) k | v G V}. 

Let furthermore 0' be the extension of <j) to G' by adding x H )■ io, y i- >• jo for some io £ I and jo G 
Consider a configuration <r : V — >• [m] with <p' C cr. For some appropriate if (cr) we have 

At(u),<t(u) = ^^ ( ^ Ar(u),<r(ii)- 
uv£E' uv£E 

Further, if cr(V') C /U J then if (<r) = A;| V|5 and otherwise, by the definition of <5 and the cell structure 
of A, we have K(a) > k\V\S + kS. In particular, as deg (Huvs-E At(u),o-(») < A • < A; we have 

The degree of cr is strictly less than fc| V|5 + A; iff <t(V) C J U J. 

We can thus compute Z^/ (</>, G) using a COUNT pin (A) oracle, by determining the values w-Na(G', 4>', w) 
for all w of degree at most k\V\5 + k - 1. This yields a reduction EVAL pin (A') < EVAL pin (A) by the 
polynomial time equivalence of COUNT pin (A) and EVAL pin (A). 

Let A" be the matrix obtained from A' by substituting X with 2. Trivially EVAL(A") < EVAL pin (yl / ) 
and EVAL(A") is #P-hard by Lemma[8l] □ 



Proof (of the Two-l-Cell Lemma). Let A G Z[X] mxm be a positive symmetric matrix containing at 



least two 1-cells. That is, A satisfies condition (A). Application of Lemmas 8.2 8.3 and 8.5 in this order 
yields the result. n 



9 The Single 1-Cell Lemma 

Lemma 9.1 (Symmetrized 1-Row Filter Lemma). Let A G Z[X] mxm be a matrix satisfying condi- 
tions (A)-(C) and containing exactly one 1-cell. Let C be obtained from A by removing all non-l-rows. 
Then 

EVALP in (CC T ) < EVAL pin (A). 
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Proof. As A satisfies conditions (A) - (C) and contains a single 1-cell we see that there is an r < m 
such that with / = [r] the principal sub matrix Ajj forms this single 1-cell. Recall that EVAL pin (^4 2 ) < 



EVAL pin (^4) by 2-stretching (cf. Lemma 5.1 ), and we have 



{A )ij — ^ AjfcAjk — ^ AjkAjk + ^ AikAjk. 

k=l k=l k=r+l 

Therefore, for ij € [r] we have (A 2 )ij = r + Y^k=r+i AikAjk, i.e. a polynomial with a constant term. 
On the other hand, if either i ^ I or j ' ^ I then A 2 - is divisible by X. Further we have C = j4[ r ][ m ] and 
thus CC T = (^4 2 )[ r ][ r ]. Hence CC T is the submatrix of A 2 which consists of exactly those entries not 
divisible by X. Recall that for any graph G and pinning <fi, we have 



and by the above 



Z A 2{^G)= Yl w-N A 2(G,cP,w) 

w&V A 2 (G) 

Z GCT {(j),G)= J2 w-N A 2(G^,w) 

wew A 2 (G) 

X does not divide w 



Filtering the weights w appropriately by means of the EVAL pm (yl) = COUNT pin (^4) correspondence, 
we can devise a reduction witnessing EVAL pin (CC T ) < EVAL pin (yl). n 



Single 1-Cell Conditions (SIC). We define some further conditions for matrices A € I J [X] rnxrn 
satisfying conditions (A)-(C). Define 8y = deg(Aij) for all i,j G [m] and let r := min{z G [m] \ 
An > 1}- That is, r is the smallest index such that 5\ r = 5 r \ is greater than zero and it exists because 
A has rank at least 2. 

(S1C-A) A has exactly one 1-cell. 

(S1C-B) A\^ = . . . = A r ~i ; *, i.e. the first r — 1 rows of A are identical. 

(S1C-C) For i G [r, m] we have 5\k < 5ik for all k G [m]. 

Lemma 9.2. Let A G Z[X] mxm be a matrix satisfying conditions (A)-(C) and (S1C-A). Then there is 
a matrix A' satisfying (A)-(C) and all (SIC) conditions such that 

EVAL p/ V) < EVAL p/n (^). 

Proof. The matrix A already satisfies (A)-(C) and (S1C-A). Observe first that if A satisfies (S1C-B) as 
well, a matrix A 1 whose existence we want to prove, can be defined as follows. For t G N, let A' = A'{t) 
be given by 

A'ij = AijiAaAjj) 1 for all i, j G [m]. 

Clearly, for all i G N the matrix A'(t) satisfies conditions (A)-(C), (S1C-A) and (S1C-B). 

We claim that there is a t G N such that A'(t) also satisfies condition (S1C-C). To see this, note first 
that the degrees of A satisfy mm{5 rr , . . . , 5 mm } > as all 1-entries of A are contained in the single 
1-cell >lr r _ 1 ]r r _ 1 i. Therefore, there is a t such that 5\j < t ■ min{<5 rr , . . . , 5 mm } for all j G [m]. Fix 
such a t and note that 

deg(^4- ■) = 5ij + t- Sa + t- 5jj. 
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Thus, with 5 n = we see that for all j G [m] and i G [r, m] 

deg(^' lj ) = 5ij + i • Sjj <t-5a + t- Sjj < deg(A^). 
This proves that A' satisfies condition (S1C-C). Reducibility is given as follows 
Claim 1. For all t G N we have EVALP in (A'(t)) < EVAL pin (,4). 

Proof. Let G = (V, E), <f> be an instance of EVAL pm (A'), Let G' = (V, E') be the graph obtained from 
G by adding t many self-loops to each vertex. Then Za'{4>, G) = Za(<P, G'), witnessing the claimed 
reducibility. -I 

It remains to show how to obtain condition (S1C-B). We give the proof by induction on m. For m = 2 
this is trivial. If m > r > 2 assume that A does not satisfy (S1C-B). Define C = Ar f ._]ir m ], i.e. the 
matrix consisting of the first r — 1 rows of A. By the Symmetrized 1-Row Filter Lemma |9T] we have 
EVAL pin (CC T ) < EVAL pin (A). Further, C has rank at least 2 as C a = 1 for all i G \r - 1) but the 



rows of C are not identical and therefore CC T has rank at least 2, as well (see Lemma 5.2 ). 



Application of the General Conditioning Lemma |4 .4 1 yields EVAL pin (C / ) < EVAL pin (CC T ) for a 
k x k matrix s atisf ying (A)-(C) such that k < r — 1. If C' has at least two 1 -cells we apply the Two 
to obtain #P-hardness of EVAL pin (C"). Therefore we can chose some A' satisfying 



4.5 



1-Cell Lemma 

the conditions ofthe Lemma such that EVAL pin (A') < EVAL pin (C")- If C has only one 1-cell then the 
proof follows by the induction hypothesis, as C' has order k < r — 1. □ 



Definition of and C^. For the remainder of this section we fix A G Z[X] mxm satisfying condi- 
tions (A)-(C) and all (SIC) conditions. Furthermore, 5ij for all i, j G [m] and r are defined for A as in 
the Single 1-Cell Conditions. For k G N define the matrix A^ by 

A [ § = A^A^f' 1 for all i,j G [m\. (9.1) 

Let further, be defined by 

c [k] =aW(aW) t . (9.2) 

Lemma 9.3. Let A and be defined as above. Then 

EVAL p/n (C [fc] ) < EVAL p/n (A). 

Proof. Let G = (V, E) and a pinning be an instance of EVAL pin (CM ). Define a graph G' = (V', £') 
as follows: 

I/' := V U {z} U {u e | e G E 1 } 

£' : = {uv e ,v e v \ e = uv <E E} U {(v e z) 2k ~ 2 \ e G 

The edges u e z have multiplicity 2fe — 2. Let 0' be the extension of <j) defined by U {z h- » 1}. Let 
a 5 0' be a configuration on G' and A Then its weight equals 



uv&E' uvGE v=\ 

= TT c [k }, <y 

11 <x(m),it(o) 



uvEE 

The last equality follows directly from the definition of CM. This yields EVAL pin (C[ fc l) < EVAL pin (,4) 
as required. n 
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We need some further notation to deal with polynomials. For / G Q[X] and A G C let mult(A, /) 
denote the multiplicity of A in / if A is a root of /, and mult(A, /) = 0, otherwise. The k-th root of 
A G C is the /c-element set X 1 ^ = {fx G C | fx k = A}. Slightly abusing notation, X 1 ^ will denote any 
element from this set. For every root A of , every r < j < m and k G N we define 

m(X,j) = min{ mult(A 1/fc , cfj) \ k>l}. (9.3) 
The following Lemma is the technical core of this section. 

Lemma 9.4 (Single- 1-Cell Technical Core). Let be defined as above. The following is true for all 
j G [r, m]. 

(1) For any root X of and all k G N we have 

(la) mult(A 1 / fc ,cf 1 1 ) = mult(A,cJ?) > m(X,j). 
(lb) mult(A 1 / fe , Cf)) > m(X,j). 

(2) If row A\ * and As * are linearly dependent then for any root X of C-q cotg? all k G N we Ziave 

mulKA 1 /*, (7}?)= mult(A 1/fc ,cfj ] ) = mult(A 1/fc , c]^ 1 ). 

(3) If row Ai t * and Aj t * are linearly independent there is a root X of such that 

mu\t(X,C$)>m(X,j). 

The proof of this lemma is technically involved and quite long. We therefore show first, how to prove 
the Single-l-Cell Lemma, provided that the above holds. 

The following basic fact^] about polynomials will be used frequently in the following. 

Lemma 9.5. Let f G Q[X] and A G C. 

(1) There exists a unique (up to a scalar factor) irreducible polynomial p\ G Q[X] such that X is a 
root ofp\. If /(A) = then p\ divides f. 

(2) If mult(A, f) = s then f = p s x f for some f G Q[X] which satisfies /(A) / 0. 

(3) Iff(X) = then X 1 ^ is a root of f(X k ) and mult(A 1 / fc , f(X k )) = mult(A, f(X)). 



Proof (of the Single l-Cell LemmaWB. Let A G Z[X] mxm satisfy conditions (A)-(C) and (S1C-A). 



By Lemma |9^2| we may assume w.l.o.g. that A indeed satisfies all (SIC) conditions. By the Two-l-Cell 



Lemma 4.5 it suffices to prove the existence of a positive symmetric matrix C which contains at least 
two 1 -cells such that 

EVAL pin (C) < EVALP in (,4). 

Recall the definition of in equation ( |9.2| ) and the definition of r in the (SIC) conditions. We prove 
a technical tool. 

Claim 1. There is a root A of cj 1 ^, an index j G [r, m] and k G N such that 



5 cf., for example, Chapter IV in [ Lan02l 
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(1) mu\t(\V k ,d$) < mult(A 1 / fe ,c{i ] ). 

(2) ForalH G [r,m] we have mult(A 1 / fc , cfj) < mult(A 1 / fc , c[ fc] ). 

Proof. Choose A and j G [r, m] such that m(A, j) is minimal for all choices of A and j which satisfy 



Al * and A^* are linearly independent and mult(A, cj 1 ^) > m(\,j). 



(9.4) 



9.4 



Note that these A and j exist by Lemma 9.4 [ 3). C hoose k G N such that, by equation ( |9.3| >, m(A, j) 
mult(A 1 / fc , c{* ] ). This proves (1) by Lemmaf" 



la). 



To prove (2) fix i G [r, m]. If and A^* are linearly independent our choice of A and j implies 
that m(X, i) > m(A, j) and therefore 



mult(A 1 / fe ,c]f 1 ) > m{X,i) > m(X,j) = mult(A i//c , C^ J ) 

where the first inequality holds by Lemma |9~4] 4b). 

If otherwise A\,* and Aj,* are linearly dependent then 

mult(A 1/fc ,c]f ] ) = mu\\{X l l k ,cf}) > m(X,j) = mult(A 1/fe , eg 1 ). 
The first equality holds by Lemma |9~4] ^2) and the inequality is true by Lemma [9~4| la). _| 

Choose j, A, k as in ClaimjlJ Define t = mill^A 1 /' , cfj). Let p\ be an irreducible polynomial such 
thatp\,(A) = 0. Let s := min{ mult(A 1 / fe , C^}) \ a, b G [m]} and define the positive symmetric matrix 



l/k fifth 



ab 

C :-- 



-C^\ 

pi 



Px ■ 



Note that s < t by definition. We consider two cases. If s = t then Cy = Cji = 1 but by Claim[TJl) 
we have C\\ ^ 1. Therefore, C contains at least two 1-cells. 

If otherwise s < t, fix witnesses a, b G [m] with mult(A 1 / fe , C^) = s. 



Claim 2. Either a > r or b > r. 



Proof. Recall that by condition (S1C-B) the first r — 1 rows of A — and hence those of A^ — are 
identical. By the definition = AM(AM) T (cf. equation dO) we have 



\2fc-2 



^=1 



^=1 



In particular, for all a', 6' G [r-1] we have C^J, = eg and hence mult(A 1 / fc , cjj,,) = mult(A 1 / fc , C*g) > 
mult(A 1 / fc ,cfJ) which proves the claim. -I 



Combining this claim with Claim 



1 ;2) and the fact that < t = mult(A 1 / fc , eg 1 ), we have either C aa + 
1 or Cm ^ 1. As C a b = Cba = 1 by the definition of a , b we see that C contains at least two 1-cells. 
We have EVALP' n (CW) < EVAL pin (A) by Lemma [93) Further EVALP' n (CM \ Px ) < EVALP' n (CW) 
by the Prime Filter Lemma[7/7l Then EVALP' n (C) < EVAL pin (A) follows from the fact that 

i G) = p s x ^Zc{4>, G) for all graphs G = (V, E) and pinnings <p. 



PA 



9. 1 Proof of the Single-l-Cell Technical Core Lemma 
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9.1 Proof of the Single-l-Cell Technical Core Lemma [94 



Proof of the Single-l-Cell Technical Core Lemma 9.4 Fix j G [r, m] and define the value b = 
mm{5ji - S u | i G [w,]}. By condition (S1C-C) we have 6 > 0. Further, if rows A^* and Aj^ are 
linearly dependent, then Aj* = X b A\^. 



For simplicity of notation, define a% := 5u, hi := <5jj — 6 and a := bi 



5 



8\i — b for all 



i G [m]. Observe that all a are non-negative and moreover a = for all i G [m] iff rows Ax,* and A^* 
are linearly dependent. If these rows are linearly independent then not all a are equal. We have 



A 
A 



L* 



■J,* 



( X ai ... X a ™ ) 
( x b+bl x b+bm ) 



By the definition of A^ in equation (9J_l, we have A 

b + b u + (k — \)a v = b + c v + ka v we have 



[k] 
iv 



A iv (A lv ) k - 1 = A iv X^- x ^ and by 



A 
A 



[k] 
1.* 
[k] 



Therefore, by = A^(A^) T , 
r [k] 

r [k] 

Cf) 



j.i 



X b ^J^ c i+2fcai _|_ 
j£2b^jr2ci+2kcti _|_ 



_|_ j^2ka m 

_|_ J^c m +2fca m 

_|_ j^2c m -\-2ka„ 



(9.5) 



Let A be a root of cf} and recall that a\ = . . . = a r -x 
(5 lr = a r > 0. Therefore A / and by equation ( |9.5| ) and Lemma 9^3) we see that \ l / k is a root of 
C [ f} with 

(9-6) 

which, by \9.5) implies 



as r was defined to be minimal such that 



mult(A 1/fe ,cf 1 ] ) = mult(A,C 

Assume first that Ax t * and Aj^ are linearly dependent. Then c\ 

C [ f) = X b ■ Cf} and cf} = X 2b ■ cf}. As A / we have 



ii / 



mu\l(\ 1/k , cf}) = mu[t(A 1/fe ,C^ 1 ) = mul^A 1 ^, cf}) 



This proves statement (2) of Lemma 9.4 and the case of Lemma 9.4 1) where A\^ and Aj t * are linearly 
dependent. 

It remains to prove statement (1) for linearly independent Ax,* and Aj^ and statement (3). Assume 
in the following that A\^ and Aj * are linearly independent. If furthermore there is a A; such that \ l l k is 
not a root of cf} then m(X, j) = and the result follows. Assume therefore that \ l / k is a root of cf} 
for all k G N. Then, for all k > 1 we have 

= Cf (X 1 ^) = \ b /k ^\2 a i+ c i/k _|_ _|_ y2a m +c m /k^ 

Defining f\(z) = A 2ai+ClZ + . . . + A 2a ™ +C "^ we get 

= Cf}(\ 1/k ) = \ b ' k f x (l/k) for all k > 1. (9.7) 
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Let denote the Z-th derivative of f\ and let a G C be such that A = e a . Then we have /a(^) = 

e a(2a 1 +c 1 z) _|_ _ _ i e a(2a m +c m z) an( j fa us 

/«(0) = (aci)'e 2aai + . . . + (ac m y e 2aa ™ = (aci)*A 2ai + . . . + (ac^A 20 ™. (9.8) 
We will prove the following property of the derivatives of f\. 

For all I > 1 we have (0) = 0. (9.9) 
The proof relies on the following claim. 

Claim 1. Suppose that g{z) = u(z) + iv{z) is a function that is analytic in the real segment [0, 1] and 
{r n } n6 N, {s n } ng N are sequences from the real segment [0, 1] such that lim n _ > . 0O r n = linin^oo s n = 
and u(r n ) = v(s n ) = for all neN. Then 

(a) g(0) = 0. 

(b) There are sequences {r' n } n ^ and {s' ra } n eN m the real segment [0, 1] such that 

lim r' n = lim s' = 

and u'(r' n ) = v'(s' n ) = for all n G N with v! , v' being the derivatives of u and v. 
Proof. W.l.o.g. we may assume that {r„} ng N and {s n } ng N are monotone. Then, since g is continuous, 

g(0) = lim g(z) = lim u(z) + % lim viz) = lim u(r n ) + i lim v(s n ). 

z— >0 z— >0 z— >0 n— >oo n— >oo 

Denote by uq, vo the restrictions of u and v to the real segment [0,1]. Then uq, vq are continuous and 
differentiable, and as u(r n ) = v(s n ) = for all n G N we see thag there are G [r n+ i,r n ] and 
s' n G [sn+i, «n] such that Uo(rjj) = v' (s' n ) = for all n G N. H 



To prove (9^9 1, recall that by equation ( |9.7| ) and the assumption that A/Owe have f\(l/k) = for all 
k G N. The definition of /a implies = f\(l/k) = u(l/k) + iv(l/k) so that the conditions of Claim[T] 
are satisfied. Applying Claim [T] inductively on the derivatives of f\ then yields ( |9.9[ ). 

In the following, it will be convenient to partition [m] into equivalence classes Nq, . . . ,N t such that 
i,j S N v for some ^ G [0, t] iff q = Cj. For each let c v be a representative of the q values pertaining 
to N u . Since not all values a are equal, we have t > 1. Further, condition (S1C-B) implies that 
ci = . . . = c r _i and we assume w.l.o.g. that Ni 3 [r — 1]. By definition, there is a = and we 
assume that A*o is its corresponding equivalence class. Define, for each v G [0, t], a polynomial 

g v {X) = x2ai - (9-10) 

ieiV„ 



By equations ( |9.8| > and ( |9.9| > we obtain the following system of linear equations, for I = 1, . . . , t 

= /l°(0) = (aco) l go(X) + ... + (ac t ) l g t (X) 
= {ac 1 ) l g 1 (\) + ... + (ac t ) l g t {\) 



Recall that this is the Mean Value Theorem 



9.1 Proof of the Single-l-Cell Technical Core Lemma 
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The second equality follows from cq = 0. The values c\ , . . . , ct are pairwise different and non-zero by 
definition. Therefore = ffi (0) for I = 1, . . . , t forms a homogeneous system of linear equations with 
a Vandermonde determinant. Hence this system is non-singular, implying <?i(A) = . . . = gt(X) = 0. As 
A is a root of cfj and C^ 1 (X) = g (X) + . . . + g t (X) we further infer g (X) = 0. 

Note that our considerations so far are independent of the specific root A of cjj . Therefore, we see 



that every irreducible polynomial g which divides C-q divides all the gi as well (cf . Lemma 



9.5 



1)). Let 



hi, . . . , h z be the different irreducible divisors of C-q such that w.l.o.g. the leading coefficients of these 
hi are positive. For each i G [z] let nii be maximal such that h^ 1 * (X) divides g v (X) for all v G [0, £]. 
Defining = /i™ 1 (X) • • • h™* (X) we see that h{X) divides each g„(X). Thus, for every v G [0, t], 
there are polynomials f v {X) = h(X)~ 1 ■ g v {X). Hence 

Cff = h(X)(f (X) + ... + f t (X)) and g u (X) = h(X) • f„(X) for all i/ G [0,t]. (9.11) 

Observe that the degree of at least one polynomial f v is positive. To see this, recall that by N\ ~D [r — 1] 
and ai = ... = a r _i = the polynomial g\ has a non-zero constant term. Furthermore, by the fact 
that all 1 -entries of A are contained in 1 -cells and A contains exactly one such 1-cell (by condition 
(S1C-A)), for all v ^ 1 the polynomial g v does not have a constant term. Therefore g\ and g v differ by 
more than a constant factor, for all v 7^ 1. With h(X) being the greatest common divisor of the g v the 
existence of such an f v now follows. The leading coefficients of all f v are positive, since that of h(X) 
is and all coefficients of the g v are positive. Therefore, 

deg(/ (X) + ... + /t(X)) >0. (9.12) 

Claim!. m(X,j) = mult(A,/i). 

Proof. It follows from the definition of m(A, j) in equation ( |9.3| ) that we have to prove that mult(A 1 / fc ,C^ ] 

mult(A,/i) for all k and that there is a A; such that equality holds. Clearly, mult(A 1/,fc , ) > mult(A,/i) 
holds for all k > 1 because, by equation ( |9.5| > and the definition of the g v in ( |9.10[ ), we have 



Cg 1 = X b (x £ °g (X k ) + ... + X £ tg t (X k ] 

= X b h{X k ) (x £o f (X k ) + ... + X £t f t {X k ] 



It remains to find a A; such that mult(A 1 / /c , cfjJ) = mult(A,/i). Define a polynomial f(X, z) = 
X^foiX) + . . . + X £tZ f t (X), then 

Cg 1 = X b /i(X fe ) • f(X k , k- 1 ) (9.13) 

Note that there is an I < t such that f(X,^) / 0. To see this, assume the contrary and note that 

/ (a, = (A £ °/' ! y/o(A) + . . . + (\ £t/tl ) l ftW 

Since the values cq, . . . , q are pairwise different, this gives rise to an invertible system of linear equations 
with a Vandermonde determinant: 

= (A £ ^ ! )/o(A) + ... + (X^ t[ )f t (\) 
= (A £ o/* ! ) 2 /o(A) + ... + (A^ ! ) 2 /*(A) 



> 



= (A £ °/* ! )Vo(A) + ... + (A e */* ! )*/*(A) 
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Therefore, /i(A) = . . . = /t(A) = which implies that all /„ have a common root. By Lemma 9.5 1) 
they thus have a common non-trivial factor, in contradiction to the definition of h. 
Hence we may fix an I < t such that /(A, ^) / 0. Define k = tl/l, then 

+ /(A, ^) = /(A, A;' 1 ) = /((A 1 / fe ) fe , A; -1 ). 

That is, A x / fc is not a root of f(X k , k' 1 ) which implies mult(A 1 / fc , cf}) = mult(A, h) by equation 



(9.13). 



Let us now prove statement 9.4 lb). By equation ( |9.5[ ), we have 

Cf} = X 2b [x 2 ^g Q {X k ) + ... + X 2 ^g t (X k ) 

= X 2b h(X k ) (x 2£o f {X k ) + ...+ X 2£t f t (X k ) 



Therefore, mult(A 1 / fc , cfj) > mult(A 1 / fc , h(X k )) = mult(A, h(X)), which, by Claim Q impl 
mult(A 1 / fe , df ) > m(A, j) — just as claimed in (lb). 



To prove (la) and (3) recall that, by equation ( |9. 1 1[ ), 

eg = h^{X) ■ ■ ■ h™°{X) (/ (X) + . . . + f t (X)) 

which implies that every root A of fo{X) + . . . + ft(X) is a root of c{^. And since the hi are the 
irreducible divisors of such a A is therefore a root of one of the hi as well. By equation ( 9.12 i and 
the fact that all roots of are non-zero, the polynomial fo(X) + . . . + ft(X) has at least one non-zero 
root, and hence 

mult^cji 1 ) > mult(A,/i). 

By equation thus mult(A 1 / fc , cf}) = mult(A, eft) > mult(A, /i) and (3) then follows by 
Claim [2] This also proves (la) for the case that this root A of C-JJ is a root of fo(X) + . . . + ft(X). 
Otherwise mult(A, C^) = mult(A, /i) which proves (la) by Claim [2] This finishes the proof of the 
lemma. 
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